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Introduction 


This document is the final report for contract 
NAS8-31771. The contract was awarded to the University of 
Arkansas at Fayetteville by George C. Marshall Space Flight 
Center, Marshall Space Flight Center, Alabama. The prime 
objective of the contract was to perform data processing 
and graphical display of Laser-Doppler Dust Devil experi- 
mental data collected by MSFC in Arizona in August, 1975. 

This report assumes familiarity with NASA Report 
No. FAA-RD-74-213 , Development of a Laser Doppler System 
for the Detection, Tracking, and Measurement of Aircraft 
Wake Vortices. The NASA report covers the theory of 
Doppler systems, signal processing, and data processing 
display. 

A summary of the research activities conducted for 
the Dust Devil contract is given in Figure 1. The remainder 
of this report contains a description of the software that 
was developed to process and display the Dust Devil data. 

The program listings are included in the appendix of this 
volume; the graphical displays are included in a separate 
volume . 

Software Development 

The University of Arkansas research team selected 
an interactive approach to analyze and display the dust 
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devil data. One of the prime reasons for selection of this 
approach is the programs are intended to be used as produc- 
tion programs and executed many times. The software was 
not developed for a single run using the batch mode of 
operation. 

Since the developed software is a research tool, 
the interactive nature of these programs allows the user to 
make qualitative decisions in a real-time mode of operation. 
Thus, the turnaround overhead inherent in batch operations 
is avoided by the interactive approach. In addition, since 
the analysis is repetitive and essentially open-ended, the 
developed software allows the user access to all the dust 
devil data. Hence, it is not necessary to resort to the 
use of "representative" or "average" data, although the user 
may do so if he desires. 

The developed software is of a modular format 
which is common to interactive systems and allows for 
expansion and refinement of the system. The system was 
designed so such modification to the system can be accom- 
plished with minimum effort. 

Dust Devil Data 

MSFC provided the University of Arkansas research 
team with fifty-eight 800-bpi tapes containing the dust 
devil data to be displayed and analyzed. These tapes contain 
approximately 200 million bytes of information collected 
using the laser doppler system. 
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These data were edited; defective data, bad tapes, 
and duplications were eliminated from the data set. Thirty- 
seven runs of the dust devil data were found usable, 
representing 559 full scans which include the minimum- 
maximum angle or the maximum-minimum angle. These 559 
full scans include 261,178 frames. Table 1 summarizes the 
usable data. 

Conversion Software 

Each dust devil data tape contains a header record 
followed by a frame data record for each frame. The header 
record consists of the run scan characteristics and the 
starting clock time of the run. The frame data record con- 
tains the frame number, the S, Y positional information in 
Cartesian coordinates, R the radius of scan field, and T 
the scan angle. In addition, the clock time and 104 values 
representing the amplitudes of the return signal for the 104 
velocity cells are recorded. 

Because the word length on these machines is 32 bits, 
and the word length on the dust devil data tapes is 36 bits, 
the data as provided were not directly acceptable to the 
IBM 360 and IBM 370 digital computers. This conversion 
difficulty was overcome by the development on the University 
of Arkansas 370/155 of an assembly language usable on other 
IBM 360-370 machines. 
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TABLE 1 



SUMMARY OF 

USABLE DUST 

DEVIL DATA 


Day 

Run 

Section 

Scans 

Frames 

234 

9 

1 

16 

4970 

234 

14 

1 

21 

6539 

234 

14 

2 

5 

1562 

234 

17 

1 

3 

951 

235 

4 

1 

35 

10879 

235 

5 

1 

25 

7774 

235 

5 

2 

4 

1246 

235 

5 

3 

3 

928 

235 

11 

1 

35 

10875 

235 

23 

1 

31 

9629 

236 

6 

1 

37 

11501 

236 

6 

2 

9 

2782 

236 

6 

3 

4 

1252 

236 

7 

1 

53 

16525 

236 

17 

1 

30 

18685 

236 

17 

2 

4 

2479 

236 

17 

3 

5 

3113 

236 

17 

4 

7 

4358 

237 

6 

1 

3 

1871 

237 

6 

2 

2 

1229 

237 

6 

3 

1 

604 

237 

6 

4 

4 

2501 

237 

6 

5 

1 

620 

237 

9 

1 

9 

5577 

237 

12 

1 

27 

16812 

238 

2 

1 

23 

14336 

238 

3 

1 

29 

18050 

238 

4 

1 

6 

3722 

238 

4 

2 

23 

14327 
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TABLE 1 continued 


Day 

Run 

Section 

Scans 

Frames 

238 

5 

1 

1 

1380 

238 

6 

1 

9 

5600 

238 

6 

2 

21 

13056 

238 

7 

1 

17 

10572 

238 

8 

1 

12 

7453 

238 

9 

1 

16 

9947 

238 

10 

1 

17 

10603 

238 

11 

1 

11 

6870 
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When called for the first time, this routine reads 
the header record and stores the fraction bits of the initial 
clock time for later use. The first call and all succeeding 
calls then reads a frame record. The integer items (frame 
number and the 104 amplitudes) are easily converted to 36 
bits by eliminating the high-order 4 bits which are always 
zero for these data. The fraction bits of the saved initial 
clock time are subtracted from the fraction bits of the 
frame clock time to yield an integral number of clockticks 
which is returned to the calling program. This leaves the 
4 floating-point items X, Y, R, T, which are handled 
identically. 

The 36-bit floating-point representation uses one 
sign bit, 8 exponent bits, and 27 fraction bits. The 
exponent is kept in excess 128 notation and indicates the 
number of bits the assumed radix point must be moved to 
obtain the actual floating-point number. The 32-bit 
floating-point representation uses one sign bit, 7 exponent 
bits, and 24 fraction bits. The exponent is carried in 
excess 64 notation and indicates the number of hexadecimal 
digits (4 bit groups) the assumed radix point must be moved 
to obtain the actual floating-point number. 

Since all quantities are positive, the sign bit is 
always zero and presents no problems. The input exponent 
bits, minus 128, are first increased if necessary to the 
next higher multiple of 4, and after adding 64, become the 
output exponent bits. In order to maintain the correct 
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radix point position, the input exponent bits are shifted 
right if necessary by the amount of the increase. The 24 
left-most input fraction bits, perhaps modified by the right 
shift, are then joined to the output exponent bits. The 
25th input fraction bit is then tested, and if it is a one 
bit, one is added to the output fraction bits. This 
accomplishes rounding and completes the conversion for 
floating-point numbers. 

Active Area Locator Program 

The volume of the dust devil data is quite large, 
and the probability that data from a given frame is useful 
for analysis and display purposes is remote. These two 
conditions precipitated the need for software that can yield 
information concerning the frames that are most likely to 
yield desired data. Such a software package was developed 
and is referred to as the active area locator program. 

The fundamental assumption underlying the active area 
locator program is that the passage of a dust devil through 
a given area represents a disturbance to the prevailing 
environmental velocity distribution for that area. This 
disturbance will be characterized by relatively high ampli- 
tude signals from velocity cells which would otherwise 
signal weakly or not at all. That is, over time there should 
be a greater than normal variation in signal amplitudes for 
areas that experience dust devil passage. The velocity cells 
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experiencing this higher variation should indicate the 
characteristic velocities for the dust devil passage. 

After initialization, input cards describing the 
scans for a run are read, and the permanent display screen 
coordinates for the location grid used by the active area 
display are calculated. The first display identifies the 
run and allows the user to set the size of the working grid 
which will be imposed upon the laser scan field. The grid 
width is expressed in whole degrees (1° - 9°) , and the grid 
length is expressed in whole meters (1 - 99) . The user 
also sets the percentage of observed variations to be 
ignored during processing (1% - 99%). This determines the 
selectivity of the locator program, and also influences the 
running time. 

The next phase of processing involves the use of the 
conversion program to read from either one or two raw data 
tapes those frames belonging to the scans selected by the 
input cards. A composite grid location value is calculated 
using the radius and angle values from the frame and the 
grid width and length values given by the user. Each frame 
lying within a given grid is assigned the coordinates of the 
near left-hand corner of that grid as viewed from the van. 

An output routine is used to write each frame with, the grid 
value added, onto a disc data set. This output data set 
is sorted using the calculated grid value as the primary 
sort field, and the frame number as the secondary sort 
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field. The resultant sorted data set consists of frames in 
chronological order for each grid in the laser scanning 
field. 

The next phase of processing involves reading the 
sorted data set and calculating for each grid value the 
observed amplitude variation for each of the 104 velocity 
cells. A frequency count is maintained for each velocity cell 
which gives the number of times a given amplitude variation 
was observed. After the sorted data set has been read, the 
amplitude variation corresponding to the user-selected 
cutoff point is determined for each velocity cell. The sorted 
data set is then read again, and an array is prepared which 
gives the velocity cell and the grid value for each instance 
of a cutoff point being exceeded. An array containing the 
total amount by which the cutoff point is exceeded is also 
prepared for each velocity cell. 

A display is then presented which has the velocity 
ells along the horizontal axis and the total amount by which 
a cutoff point for a velocity cell is exceeded along the 
vertical axis. This display indicates the velocities 
which can be considered most active in the sense of larger 
variations in signal amplitudes over time. A hard copy of 
this display is usually made for future reference. 

The user then has the option of: (1) terminating 

the program, (2) selecting another velocity cell, (3) chang- 
ing the grid size and/or rejection percentage, or (4) pro- 
ceeding to the 3-dimensional display routines. If option 
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(4) is chosen, the user defines the boundaries of the region 
to be displayed (minimum and maximum radius and angle) . 
Control is then passed to the 3-dimensional display routines. 
When control is returned, the user is again presented with 
these four options. The subroutines used by the active 
area locator program are provided in the appendix. 

Three-Dimensional Display Programs 

The 3-dimensional software is interfaced with the 
active area locator program by subroutine MENU. The 
programs are invoked only after a region containing sig- 
nificant amplitude variations is identified in the laser 
scanning field. This technique eliminates the detailed 
analysis of insignificant data and increases the efficiency 
of the analysis procedure. 

Subroutine MENU offers three options to the user: 

(1) return to the active area locator program, (2) process 
a full scan of data with subroutine DISPLY, or (3) process 
a particular angular region with subroutine RAY. After 
return from DISPLY or RAY, subroutine MENU cycles, again 
offering options (1), (2), and (3) to the user. 

Subroutine DISPLY operates at the scan level. The 
first screen display presents the user with several 
options. Other than the option allowing for return to 
subroutine MENU, these options determine the quantity to 
be associated with the vertical axis in a 3-dimensional 
display. The base quantities are always radius and angle. 
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thus displaying the selected quantity as a function of 
position. The selectable quantities are: (1) maximum 

amplitude across all velocity cells, (2) amplitude at the 
velocity cell chosen in the active area locator program, 

(3) the velocity cell associated with maximum amplitude, 

(4) maximum velocity cell with amplitude above a threshold, 
and (5) the product of (1) and (3) . 

If option (4) is selected, the user may enter an 
amplitude threshold. For all display options, the required 
quantities are calculated by processing the frames 
associated with a given scan. A uniform position grid is 
imposed upon the laser scanning field, and the desired 
3-dimensional graph is plotted by calling subroutine 
THREED. 

Subroutine DISPLY then cycles, allowing the user to 
choose various options before returning to subroutine MENU. 

An additional option available after the first pass through 
DISPLY is a call to subroutine RAY. It is unnecessary to 
reprocess the data frames since all quantities are calculated 
during the first pass. 

Subroutine RAY is used to display a graph of ampli- 
tude over a base of velocity and radius, with the input 
data coming from a particular ray of the laser scan. As 
such, this is a very detailed analysis and is normally 
undertaken only for angles which show promise of exhibiting 
dust devil passage. 
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Subroutine RAY is cyclic in nature and works within 
the angular bounds established by the active area locator 
program. The user has the option of returning to the 
calling program, processing the entire angular region, or 
specifying a smaller angular region which includes a single 
angle as a special case. Before each ray display is pre- 
sented, the user is given the option of bypassing the 3- 
dimensional graph. After a given angular region has been 
processed, the RAY subroutine cycles, allowing the user to 
select another range of angles. The THREED subroutine is 
used for the 3-dimensional display itself. 

Subroutine THREED is an interface routine connecting 
the display routines to the 3-dimensional surface plotting 
package developed by Dr. D. L. Beveridge at City University 
of New York. The primary function of THREED is to establish 
the viewpoint coordinates for displaying the surface, and 
to define the other constants required by the CUNy package. 
The subroutines used by the 3-dimensional display programs 
are lised in the appendix of this report. 
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The following pages contain Tektronix hardcopies taken 
during a session with the Active Area Locator Program. This 
session was followed immediately by a session with the 
3-Dimensional Display Programs. The hardcopies for that 
latter session are in the back section of this appendix. 

Hardcopies 1, 2, and 3 show the input parameters 
being supplied to the Active Area Locator Program in the 
order of grid width, grid length, and cutoff percentage, 
and Hardcopy 4 shows results from the first phase of process- 
ing. The velocity cells lie along the horizontal axis, and 
the vertical axis indicates the amount by which the 90% 
point was exceeded by the observations surviving the cut. 

Hardcopies 5, 6, and 7 indicate a change in the cutoff 
point from 90% to 95%, -while Hardcopies 8 and 9 show a 
change in grid dimensions to 1 degree by 2 meters. Hardcopy 

10 shows a request to examine velocity cell 61, with Hardcopy 

11 being the result. This display indicates the most 
extreme 5% of the observations for velocity cell 61, by 
location, on a polar grid. 

Hardcopy 12 indicates how the user requests another 
velocity cell display, and Hardcopy 13, 14, and 15 are 
displays for velocity cells 62, 47, and 50 respectively. 

Hardcopy 16 indicates a request to change input parameters, 
with Hardcopy 17 showing the cutoff point increased to 98%. 
Hardcopies 18 through 22 show the differences in the displays 

resulting from this input parameter change. 

Hardcopies 23 and 24 are preparatory to entering the 
3-Dimensional Display Programs and indicate the region to 
be analyzed. 
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ACTIVE AREA LOCATOR PROGRAM 

DIMENSION I AMP (104 ), I MIN (104)/ I MAX (104 1, I FREQ ( 2 50, 104 ) 

DIMENSION 199(104), IC OUNT ( 1C 4 ) , 10 V ER ( 10 4 ) > POS ( 4 ) , I VAR (10 4) 
DIMENSION I DAT A (20000, 2 >, LINE 1(41), LI NE2(38)»L I NE 3 ( 6 ) * L I NE 4 ( 58) 
DIMENSION JKEE PK20CC0) 

DIMENSION L INE 5( 21 ) ,L INE6( 48 ),L INE7( 76) 

DIMENSION I START (60 ), T 1 ( 60 ), I STOP ( 6C ) , T2 ( 60 ) 

DIMENSION ITITLE( 36),L INEb(46),LINE9(47),LlNE13(64), ITAPES(2) 
DIMENSION L INE 11 (45 ) 

DIMENSION ISP0TX(U,i9),ISP0TY(i:,19) 

DIMENSION MINI (3), MAXI (3), l INE 12(32 >»L INE 13 (50) 

DIMENSION L INE14C6) 

EQUIVALENCE ( X,PCS( 1 ) ), ( Y,PGS(2) >, (R,PGS(3 ) I, (T,POS(* I) 

DATA LINE1/86, 69, 76*7 9,67, 73, 84, 89, 32,67, 69, 76, 76,83, 32, 65, 78,68, 
132,65,77,7 9,65,78,8 4, 8 3, 32,79,86,6 9,8 2, 32, 32,32,37, 32,80,79, 73, 7 8, 


184/ 

DATA LINE2/82, 69, 67, 79, 82, 68, 32, 67, 69, 76, 76, 83, 32, 7 9, 70,32, 73,78, 
184,69,82,69,83,84, 32, 65, 76 , 6 8, 32 , 69, 78 , 84 , 69, 82, 22, 3 4, 68, 34 / 

DATA L INE 4/ 69, 78,84,69, 82,32,51,32,68,73, 71,73, 84, 32, 67, 69, 76,76, 
13 2, 78,85,77,66,6 9,8 2, 32, 79, 8 2,32,34,83,34,32,84,79, 32,83,84,79, 

18 0, 32, 79,82, 32,3 4,71, 34, 32, 84,79, 32, 67, 89, 6 7, 7b, 69, 5 8, 32 / 

DATA LINE 5/ 69, 78, 84,69,82, 32,34,68,34,32, 8 4,79, 32,6 7, 79, 78,84,73# 
178, 85,69/ 

DATA LINE6/69, 8b, 84, 82 , 6 9, 77 , 69, 32, 80, 79, 73,78,84,83,32, 70,79,82, 


132,86,69,76,79,67, 73, 84 , 69, 32, 67,69,76, 76, 32/32,32,32,32,32,32,32, 
132, 80,69,82, 67 ,69, 78, 84 / 

DATA LINE 7/69, 78, 84, 69, 82, 32, 34, 83, 34, 52, 84, 79, 32,8 3,84,79,60,32, 
134, 86, 34, 32, 70,79, 82, 32, 65, 32, 86, 69, 76, 79, 67, 73, 84, 89, 32,67, 69,76, 
176, 32,34,71, 34,3 2, 84, 79, 32,6 7, 72, 6 5, 7 8, 71,69, 32, 71, 82,73,63,32,7 9, 
18 2, 32,34,77, 34,32, 70, 79, 82,32,77,6 9,78,85/ 

DATA LINE10 /69,78, 84, 6 9, 82, 32 » 80, 69, 82 , 67 , 69, 78 , 8 4, 32, 79, 70, 32,7 9, 
166,83,6 9,82,36,6 9,6 8, 32, 86,65, 82, 73,65, 84 , 73, 7 9, 78, 8 3 , 32, 34, 79, 32 , 
166,69,32,73,71,78,79,82,69,68,44, 32,84, 87,79,32,68,73,71,73,84,83, 
158,32/ 

DATA LINES /b9, 78, 84, 69, 8 2, 32 , 71 , 02 , 73 , 6 8, 32,87, 73, 68, 84, 72, 32, 73, 
178,32,8 7,72, 79,76,6 9, 32 , 68, 6 9, 7 1 , 82, 69, 69, 83,4 4,32, 79,78,6 9,32,6 8, 
173,71,73,84,58,32/ 

DATA LINE9/69, 78, 3h, 6 9, 82, 32,71,82, 73,6 8, 32 , 76, 69, 7 8, 71, 84, 72, 32 , 
173,78,32,87,72,7 9,76,6 9,32,77,69,84,69,82,8 3,44,32, 84,87,79, 32,68, 
173,71,73,34, 83,58,32/ 

DATA IT I TLE/ 68,65, 89, 32, 32, 32,32, 32,82, 85,78,32,32, 32,32,34,65,80, 
169, 32,32,32,32,32, 32,32,65,78,68,32,32,32,32,32,32, 32/ 

DATA L INE 11/ 7 1,82, 73, 68, 3 2, 67, 7 3, 6 8, 84 , 72, 32,32,32, 32, 68,69,71,32, 
169, 69, 83, 32, 32,7 1, 82, 73, 68, 32, 76, 69,78, 71 , 84, 72, 32, 32,32, 32, 77,6 9, 
184,69,82,83, 32 / 

DATA LINE 12/ 69, 76, 84, 6 9, 82, 32, 32, 32, 32, 73, 77, 85, 77, 32 , 82 , 65, 68, 73 , 
18 5,83,32,73,78,32,87,72, 79,76,69,32, 77,69, 84,69,82,83,44,32,84,72, 
182,69,69,32,68,73, 71,73, 84,83, 58,32/ 

DATA L INE 13/ 69, 7 8, 84, 69,82, 32,32, 32, 32, 73, 77,85,77,32,65, 78, 71,76, 
16 9, 32* 73,7 8,32,8 7,72,79, 76,6 9,32,6 8,6 9,71,82,6 9,69, 83,44,32,84,8 7, 
179, 32,63,73,71 ,73,84,83, 58,32/ 

DATA L INE14/ 83, 64, 79, 82, 65, 71,69,32,69, 88,67,6 9, 6 9,63,69,68/ 

DATA MINI/77,73, 78/ 

DATA MAXI/77,65, 88/ 
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C INITIALIZE FOR CARO READING 0F QUALITY 

CALL INITTC120) 

CALL TERM( 3*1024> 

CALL C HRS I Z ( 3 ) 

I TAPE-0 
I TAPE 5 ( 2 ) *Q 
NTAPE-C 
I SCAN-1 

C READ CARO WITH SCAN LIMITS DAY RUN AND TAPE NUMBER 

54 READ! 5* 101, END *53) I S TART ( I SCAN >* Tl( ISC AN I * ISTOPI ISCAN >* T2f I SCAN >, 
1 IDAY* IRUN* I TAP El 

101 FORMAT! 1 6* F 6 . 3 * I 6* F 6 . 3* 13* 12* 15) 

C COUNT THE NUMBER UF SCANS READ 

I SCAN- 1 SCAN «■! 

C SAME TAPE AS BEFORE? YES TO READ AGAIN 

IF! ITAPE.EO. ITAPE1) GO TO 54 
C NO* COUNT ONE TAPE 

NTAPE-NTAPE *1 

C SAVE FIRST SCAN NUMBER FOR THIS NEW TAPE 

I SC AN1 * ISC AN-I 

C PUT TAPE NUMBER INTO STORAGE 

ITAPESINTAPE )* IT APE 1 

C NOW MAKS NEW TAPE LOOK LIKE OLD ONE 

ITAPE-ITAPEi 
C CONTINUE READING 

GO TO 54 

C RECORO TOTAL NUMBER OF SCANS 

53 ISCAN- ISCAN-i 

C DUMMY VALUE TO BE SURE LAST SCAN IS PROCESSEO CORRECTLY 

ISTART1 ISCAN + i 1*0 

C PLACE DESCRIPTIVE INFO INTO TITLE 

I IDAY- IDAY 

CALL DIGITS! II DAY, ITITLEI 5)>3> 

I I RUN* I RUN 

CALL DIGITS ( II RUN? ITITLE!13)*2) 

I ITAPE-ITAPES! 1) 

CALL DIGITS ( I I TAPE* I TITLE(21 )* 5) 

I ITAPE*ITAPES( 2) 

CALL DIGITS ( 1 1 TAPE* I T I T LE ( 31 ) * 5 ) 

C SET LENGTH OF TITLE DEPENDING ON NUMBER OF TAPES (1 OR 2) 

LENGTH-20 

IF! NT APE. EQ. 2) LENGTH-36 

LOOP TO 43 CALCULATES INTERSECTION OF EVERY 5 DEGREES WITH 
EVERY 50 METERS FROM 0 TO 9G DEGREES AND 0 TO 500 METERS 
SAVE THE SCREEN COORDINATES IN ISPOTX AND ISPOTY 
DO 48 1-1*11 
RR*50* ( I-l ) 

00 48 J * 1 * 1 9 
T T-5* ( J-i ) *. Q1 74533 
TX-RR*SIN!TT) 

TY»RR*COS( TT) 

TX-TX*1 . 5+40 
T Y* T Y*1 . 5+2G 
ISPOTX! I* J)*TX 
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48 I SPOT Y ( I # J ) *TY 
C CLEAR SOME COUNTING ARRAYS 

65 00 1 1*1*104 
I OVER ( I > *0 
ICOUNT ( I ) *0 

00 1 J*l* 256 
1 I FREQ ! J* I ) *0 

C DISPLAY TITLE ON NEW PAGE 

CALL NEWPAG 

CALL ANSTRCLENCTH, ITITLE > 

CALL NEWLIN 

C ASK FOR GRID WIDTH 

CALL ANSTRC46* LINES > 

55 CALL BELL 
CALL TINSTR < 1* IWIOTH) 

I WlDTH*IWlDTH-48 
IFC IWIDTH.LT. 1 .OR. IWIDTH.GT. 9) GO TO 55 

C PUT WIDTH INTO LINE11 

JWIDTH-IWIDTH 

CALL DIGITS! J W IDTH* L I N£ 1 1 C 12 ) * 2 ) 

CALL NEWLIN 
C GET GRID LENGTH 

CALL ANSTR (47* L I NE9 1 

56 CALL BELL 

CALL TINSTR ( 2» LINE3 ) 

I LNGT H*L INE3 ( 2 )-4S*lt*C L INE3! U-48 > 

IFC ILNGTH.LT. 1 .OR. ILNGTH.GT. 99) GO TO 56 
C NUMBER OF GRIDS WITH SAME ANGLE IS 500 DIYIOEO BY LENGTH 

ISPAN-50C/ILNGTH 
C PUT LENGTH INTO LINE11 

J LNGTH* I LNGTH 

CALL OIGITS ( JL NGTH* LINE11C 36 ) > 2 ) 

CALL NEWLIN 

C GET REJECTION PERCENTAGE 

CALL ANSTRC64* LINE10) 

57 CALL BELL 

CALL TINSTR (2* LINE3 ) 

IPERCT-L INE3C2 >-48 *10*1 L INE3C 1 )-48 ) 

IFC IPERCT.LT. 1 .OR. IPERCT.GT. 99) GO TO 57 
C PUT THIS INTO LINE1 

JPERCT-IPERCT 

CALL DIGITS ( JP £RCT*L I NEi C 33 >*2 > 

C AND ALSO INTO LINE6 

JPERCT-IPERCT 

CALL DIGITS! JP£RCT*LINE6C39)*2> 

C CLEAR PAGE AND INITIALIZE FOR TAPE REAOING LOOP 

CALL NEWPAG 
IRET-0 
IHOP-O 

C LOOP TO 63 READS RAW DATA TAPE ADDS GRID NUMBER AND WRITES OUT 

DO 63 I*i» ISCAN 

C LIMITING FRAME NUMBERS FOR THIS SCAN 

1 HUNT* I STAR TC I > 

I END* ISTOPC I > 
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C IF IN T He SECOND TaPE* GO AROUND FIRST TAPE STUFF 

IFINTAPE.EQ.2. AND. I.GE.ISCAN1) GO TO 6i 
C HAVE WE ALREADY READ FIRST FRAME? YES* GO AROUND 

I F ( IHOP.ta.l ) GO TO 60 
C READ FORWARD A FRAME 

58 CALL GET< IFRAME* POS* INC * UMP* IRET) 

C ARE WE AT START OF SCAN? NO* CONTINUE LOOKING 
I F I IFRAME.NE. I HUNT I GO TO 58 
C YES* DO NOT LOSE GOOD DATA 

GO TO 59 

C READ A FRAME FROM SCAN 

60 CALL GET ( IFRAM E* POS* INC* IAMP* IRET > 

C GET INCREMENT OF RADIUS 

59 IR • R / ILNGTH 

C MAKE SURE UNEXPECTED NUMBERS DO NOT PRODUCE STRANGE THINGS 

IF(R.LT.C) IR-C 
IF(R.GE.500) I R * I S PAN-1 
IFUR.Gc.ISPAN ) IRMSPAN-i 

C GET ANGLE INCREMENT *** TAKE CARE OF DOUBLE ANGLE PHENOMENON *+* 

IT-<T*2)/IWIDTH 

C MAKE SURE BAD DATA DOES NOT SINK THE SHIP 

IFIIT.LT. 0) IT-G 
IFfIT.GT.89) IT-89 

C NOW WE HAVE THE ALL- I MPQRTANT GRID NUMBER 

IGR IU-ISPAN*IT ♦ I R + 1 

C PUT OUT FRAME INFO WITH GRID NUMBER ADDED 

CALL PUT ( I GR I D * IFRAMc*X* Y* R * T* I NC * I A MP* I RE T ) 

C ARE WE DGNE WITH SCAN? NO* CONTINUE READING FRAMES 

IF( IFRAME.NE. I END) GO TO 60 

C DECIDE WHETHER OR NOT TO READ FORWARD TO FIND NEXT SCAN START 

IHQP-1 

I Ff IEND.NE. 1ST ARTt I+i ) . OR . ( NTAP E . E Q . 2 . AND. I +1.EQ. ISCANi ) ) IHOP-O 
C AND HOP THE PROCESSING FOR SECOND TAPE 

GO TO 63 

C PROCESSING FOR TAPE TWO USES GET2 SUBROUTINE. SAME PROCESSING 

61 I F ( IHOP.EQ.l ) GO TO fc7 

68 CALL GET2( IFRAME*POS* INC* IAMP* IRET ) 

I F I IFRAME.NE. I hUNT) GO TO 68 

IFI I. ED. ISCANI .AND. ISTOPf ISCAN1-1 >.E0. IHUNT) GO TO 67 
GO TO 69 

67 CALL GET2C IFRAME*PQS* INC* IAMP* IRET ) 

69 IR-R/ILNGTH 
IF(R.LT.O) I R * G 

I F ( R . GE .500 ) I R * I S PAN- i 
IFf IR.GE.ISPAN ) IR* IS PAN-1 
IT*(T*2)/I WIDTH 
IFf IT.LT.O) IT-0 
IFfIT.GT.89) IT-89 
IGR ID- ISPAN*IT*IR+1 

CALL PUT ( IGR ID* I FRAME *X, Y*R*T* INC* IAMP* IRET) 

IFf IFRAME.NE. I END) GG TO 67 
I HOP- 1 

IFf IEND.NE. I START! I +1 ) ) IHOP-O 
63 CONTINUE 
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C SIGNAL ROUTINES TO CLOSE FILES 

IRET-i 

CALL PUT ( I GR I D # IFRAME»X#Y»R#T# INC# I AMP# I RET) 

CALL GET< IFRAME#POS# INC# IAMP# IRET) 

IF! NT APE. EQ. 2) CALL GET2 ( I FR AME# PQS# INC » I AMP# I RET ) 

ASK FOR OUTPUT FILE TO BE SORTED 
PRIMARY FIELD IS GRIG NUMBER 
SECONDARY FIELO IS FRAME NUMBER 
CALL SQRTIT 

C . PROCESSING TO STATEMENT 11 R EADS SORTED DATA 
C INITIALIZE SOME SWITCHES 

IST0P1«0 
IRET«0 

C READ IN THE FIRST RECORD 

CALL GETS! I GRID# IF RAME# X# Y #R # T> I NC #.I AMP # I RE T ) 

C MAINTAIN A COPY OF THE DATA FRAMES 

8 CALL PUT(IGRID#IFRaME#X#Y#R#T#INC#IAHP#IRET) 

C IF WE ARE TO STOP (END OF FILE) GO TO 11 

I F ( ISTOPi.EQ.l ) GO TO 11 

C OTHERWISE WE HAVE THE FIRST FRAME BELONGING TO SOME GRID 

NCASE-1 

C SAVE THIS GRID NUMBER 

I GR I D1 * I GR I D 

C ONLY OBSERVATION SO FAR IS CERTAINLY MINIMUM AND MAXIMUM 

DO 4 I »1 # 104 
IMINI I ) » I AM P ( I ) 

4 I M A X ( I )» IAMP( I ) 

C READ ANOTHER FRAME 

6 CALL GETS( IGRI D# IFRAM£#X#Y,R#T# INC# IAMP# IRET) 

C WAS IT END OF FILE INSTEAD? YES# GO TO 13 

I F ( IRET.NE.O) GG TO iO 

C HAVE WE CHANGED GRID NUMBERS? YES# GO TO 7 

IF( IGRID.NE.IGRI01) GO TO 7 
C NO# SAFE TO OUTPUT THE FRAME NOW 

CALL PUT ( IGR ID#IFRAME#X#Y»R#T# INC# I AMP# IRET ) 

C UP THE OBSERVATION COUNT FOR THIS GRID 

NCASE»NC ASE+1 

C FINO NEW MAX AND MIN FOR VELOCITY CELLS 

DO 5 I * i # 10 4 

I F ( IAMP( I >.LT. IMINt I ) ) I M I N ( I ) » I AMP ( I ) 

5 IF( IAMPI I).GT. IMAX( I ) ) I MaX ( I ) • I AMP ( I ) 

C KEEP READING 

GO TO 6 

C NOTHING TO RECORD IF ONLY ONE OBSERVATION 

7 IFINCASE.LT. 2) GO TO 8 

DO 9 I *1# 1 34 

C OTHERWISE OBTAIN THE VARIATIONS 

J-IIMAXC I) — IMINT m+1 

C AND KEEP A COUNT OF THE OBSERVED VARIATION 

I VAR < I ) * J-l 

9 IFREQI J# I )« IFREQ ( J# I M-l 

C PUT OUT THE VARIATION INFORMATION BEHINO FRAMES FOR THIS GRID 

CALL PUT ( IGRID1#0#0#C#0#C#NCASE#IVAR#0) 

C AND GO PUT OUT FRAME THAT SIGNALED GRIO CHANGE 
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GO TO 9 

C PROCESS END OF FILE AS A GRID CHANGE# 8UT 00 NOT READ AGAIN 

±0 I STOP I « i 
GO TO 7 

C FIND THE TOTAL NUMBER CF OBSERVATIONS 

C NOTE THAT ANY VARIATION GETS RECOROEO SOMEWHERE# EVEN ZERO VARIAT1 

11 I SUM*0 

DO 12 I*i#256 

12 I SUM* ISUM* IFRE QC I#1 ) 

C NOW FIND THE COUNT THAT GOES WITH USER- S E L EC T E D CUTOFF POINT 

I SUM* ( ( l SUM* l P ERC T I / 1 OC ) +1 

USE THIS TU LOCATE THE CUT POINTS FOR EACH VELOCITY CELL 
AND SAVE THEM FOR LATER 

00 14 I«1#104 
J SUM*0 

DO 13 J» 1 # 2 56 
J SUM* JSUM* IF RE Q U# I > 

1 F < JSUM.GE . ISUM) GO TO 14 

13 CONTINUE 

14 I99(I)*J 

PROCESSING TO STATEMENT 17 READS SORTED DATA AGAIN 
AND RECORDS ALL INSTANCES OF CUTOFF POINTS EXCEEOED 
INITIALIZE SWITCHES 
I RE T*0 
IKEEP-C 
C READ A FRAME 

15 CALL GETS! IGRI Cl # I FRAME# X# Y, R, T# NC AS E# IVAR# IRET ) 

C END OF FILE? YES# GO TO 17 

IF! IRET.NE.G) GO TO 17 

C IS IT A VARIATION RECORD? NO# GO READ AGAIN 

I F I IFRAME.NE.O ) GO TO 15 
C YES# LOOK AT E AC F CELL 

DO 16 1*1# 1 C 4 
J * I VAR ( I ) 

C VARIATION LESS THAN CUTOFF IS IGNORED 

IFIJ.LT. 199(1) ) GO TO 16 
C BUMP COUNTER 

IKEEP* IKEEP + 1 

C STURAGE EXCEEDED? YES# GO HANOLE IT AT 64 

IF( IKEEP. GT.2GC0C) GO TO 64 
C RECORD GRID VALUE 

I DAT A I IKEEP#!) -IGRID1 
C AND CELL NUMBER 

I DAT A ( IKEEP # 2 ) * I 
C COUNT ONE MORE OVERAGE 

ICOUNTI I > ■ I C OUNT (I )♦! 

C AND SUM THE OVERAGES 

I OVER (II* 10 VER (I)U«I99( I) 

16 CONTINUE 

C GO GET ANOTHER FRAME 

GO TU 15 

C SIGNAL FOR CLOSE FILE 

64 I RET* 1 

CALL GETSC IGRI Dl# I FRAME# X# Y# R# T# NCAS E# IVE R# IRET ) 



66 CALL ANSTR! lb# LINE14) 

GQ TO 65 

FIND GREATEST OVERAGE FOR SCALING OF DISPLAY 

17 JMAX-0 

DO 16 I = l> 1 C 4 

18 IFI IOVER(I) .GT .JM4X) J M A X« I OV ER I I ) 

PUT DESCRIPTIVE INFO ON NEW SCREEN 
CALL NEWPAG 

CALL ANSTR! 41, LINE1 > 

CALL NEWLIN 

CALL ANSTRILENGTH, ITITLE ) 

CALL M0VABS!60C,765> 

CALL ANSTR(45, LINE 11) 

CALL M0VABS(6CC#74C) 

CALL ANSTR138, LINE2) 

DISPLAY DIGITS FOR HORIZONTAL AXIS 

DO 19 I»i,19 

LINE3I2)«48 

IFI <1-1 I/2>*2) .Nt.O) L I NE3 < 2 ) * 52 
L INE 3 ( 1 > ■ ( I/2)+48 
I X* 78*45*! 1-1 ) 

CALL HOVABS ( IX, C ) 

19 CALL ANSTR ( 2, L INE3 I 

DRAW LINES VERTICALLY FOR EACH CELL REPRESENTING OVERAGE 
DO 2C I * 1/ 1 C 4 
1 X>49+9*(I-i > 

CALL MOVABSUX/15) 

IY-16*(75D*I0VER!I > )/ JMAX 

20 CALL DRWABS! IX,IY) 

DISPLAY DIGITS FOR VERTICAL AXIS AND ORAW DECTILE LINES 

DO 25 I*i,lG 

II»!<II-l)*JMAX)/iC) 

IFI I I. GT. 9999) GO TO 22 
CALL DIGITS! II,LINE3,4) 

GO TO 24 

22 DO 23 J J*l, 4 

23 LINE3!JJ)-42 

24 l Y«1 5 * < I -1 ) *75 
CALL MOVaBSCG, IY) 

CALL ANSTRI4/L INE3) 

CALL HOVABS ( 47 , I Y ) 

25 CALL DRWABS ( 97 8, IY) 

WAIT FOR USER TO SIGNAL WITH A "D" 

26 CALL BELL 

CALL TINSTRC, LINE3) 

IFILINE3ID.NE .68) GO TO 26 
NOW GET A CLEAR SCREEN 
49 CALL NEWPAG 

LET USER DECIDE WHAT TO DO 
CALL ANSTR < 58, L I NE4 ) 

28 CALL BELL 

CALL T I NSTR < 3, L I Nc3 ) 

CAN GO TO STOP 

IFILINE3I1 ).60.83) GO TO 5C 
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C CAN GO CHANGE GRID SIZE AND/OR PERCENT 

IF(LlNE3(i).EQ.7l) GO TO 65 
C OTHERWISE GET SELECTED VELOCITY 

DO 27 1-1,3 

27 L INE3 < I )-L INE3 I I )-48 

IC€LL»iOO*L INE 31 1 )*lu*L INE3I 2 ) *L INE3I 3 ) 
IF( ICELL.LE.O.OR. ICELL.GT.104) GO TO 28 
C PUT INTO LINE6 

JCELL- ICELL 

CALL OIGITSI JCELL,LINE6!34),3) 

C FIND ALL OVERAGES FOR THIS CELL 

NSTORt-O 
DO 32 I-.,IKEEP 

IF( IOATAI I,2).NE. ICELL) GO TO 32 
NSTORE-NSTORE*! 

IFINSTORE.GT.2COOO) GO TO 66 
JKEEPi ( NSTORt ) -I 
32 CONTINUE 
C PUT DESCRIPTIVE 

CALL NEWPAG 
CALL MOVABSI 60C, 765 ) 

CALL ANSTRI48, LINE6) 

CALL MOV ABS I 60C, 750 ) 

CALL ANSTRILENGTH, ITITLE ) 

CALL MQVABS I 60C, 735 ) 

CALL ANSTRI45, LINE 11) 

C DRAW THE RADIUS AND ANGLE LOCATION GRID 

DO 34 I»i,19 
IX- ISPOTXt 1, I ) 

I Y - IS POT Y 1 1 , I ) 

CALL MOVABS ( IX, I Y) 

I X- ISPOTXI 11^1 ) 

I Y-ISPOTYIU, I ) 

34 CALL DRWABS ( IX, IY) 

DO 45 1-2, li 
I X» I SPOTX ( I , 1 ) 

I Y-ISPOTYI I,i ) 

CALL MOVABSI IX, IY) 

00 45 J-2,19 
IX-ISPOTXI I, J) 

I Y-ISPOTYI I, J) 

45 CALL ORWABS I IX, I Y) 

DO 42 I»i*li 

1 I* I 1-1 ) *50 
IX-ISPOTXI I , i ) -40 
IY-ISPUTYtl,:) 

CALL MOVABS I I X , I Y ) 

CALL DIGITS! II ,L INE3, 3 ) 

*2 CALL ANSTRI 3,L INE3) 

DO 35 I -i, 1 9 
11-5*11-1) 

I X- I SPOTXI 1 i, I )*10 
I Y« ISPOTYI 1 1, I )U0 
CALL MOVABS I I X , I Y ) 


INFO ON NEW SCREEN 


original Pag b * to 

OF IS 

^OOR QUALITY 
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ORIGINAL PAGE IS 

CALL OIGITS ( II >L INE3*2 ) OF POOR QUALITY 

35 CALL ANSTR ( 2* L INE3) 

C DISPLAY THE POINTS ON THE RADIUS AND ANGLE GRID 

DO 3fc I I «1» NST ORE 
J-JKEEP1 (II) 

JGR ID * IDAT A ( J» 1 ) 

JT«CJGRID-1)/ISPAN 
JR«JGR IO-ISPAN*JT 
J R* JR* I L NGT H 
JT» JT*IWIDTH 
T JT-JT*.C17<*533 
J X* JR*CQS( T JT ) 

J Y* JR*S I N( T JT ) 

JY«( JY*3)/2 
JX*( JX*3>/2 
J Y* J Y*4C 
JX» JX+20 

CALL MOVABSl JY, JX> 

L INE3( U-46 

CALL ANSTR(1*L INE3) 

36 CONTINUE 

C WAIT FOR USER TO SIGNAL FINISHED 

CALL M0V48S( 825»i5) 

CALL ANSTR ( 21 » L l NE 5 ) 

37 CALL BELL 

CALL TINSTRd* LINE3 ) 

IFCLINE3(l).NE.6d) GO TO 37 
C LET USER DECIDE WHAT TO DO AGAIN 

33 CALL NEW PA G 

CALL ANSTR! 76, LINE7) 

21 CALL BELL 

CALL TINSTRIL# LINE3 ) 

C CAN CHOOSE TO STOP 

IF1L INE3I1 ) .EQ.63) GC TO 50 
C CAN SELECT ANOTHER VELOCITY 

IF(LINE3(1 I.EQ.86) GO TO *9 
C CAN CHANGE GRID SI2E AND/OR PERCENTAGE 

IF!L IN£3(1) .60 .71) GC TO 65 
C OR CAN GO TO 3-D ROUTINES 

IF1LINE3U >.Nc .77) GO TO 21 
C GET LIMITS FOR REGION FROM USER 

CALL NEWPAG 
LINE12(7)»MINI !1 ) 

L INE12 ( 8 )*M IN I (2 ) 

L INE1 2 1 9 ) «M IN I (3 ) 

CALL ANSTR(52> LINS12) 

CALL BELL 

CALL T I N $TR ( 3* LINE 3) 

IRMIN»1CC*!L IN £3(1 ) - 48 > ♦ 10* ! L I NE3 ( 2 J -4 8 ) *L lNE3(3)-48 
CALL NEWLIN 
LINE12I7)«MAXI 11 ) 

L INE12(8)«MAXI (2) 

L INE12 ( 9 )«MAX I (3) 

CALL ANSTR! 52* LI NE12) 
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CALL BELL 

CALL T INSTR (3# L I NE3 ) 

IRMAX»10C*(LI.'4f : 3( l ) -48 H-10* < L I NE 3 ( 2 ) -4 6 ) *L I Nt 3( 3 I -4 8 
CALL NEWLIN 
L INEi 3 ( 7 )» MI NI (1) 

L INEi 3 ( 6 )* Ml NI 121 
L INEI 3 ( 9 )■ MINI (3 ) 

CALL ANSTk(50# LINE13) 

CALL BELL 

CALL TINSTR12# LINE3) 

ITMIN»lC*<LINE3<i>-48)+L INE3( 21-48 
CALL NEWLIN 
L IN£13(7)«MAXI <1 ) 

L INE13IQ )«MAXI (2 ) 

LINEi3(9l*MAXI (3) 

CALL ANSTRI5Q# LINE13) 

CALL BELL 

CALL TINSTRI2# LINES) 

I TMAX«iO*IL INE 3( 1 )-48 )+L INE3I2 )-48 
C HAND CONTROL TO 3-D DISPLAY PACKAGE 

CALL NEWPAG 

CALL MENU ( I STA RT # Ti# ISTOP# T2# ISCAN# I SC AN1# I OA Y# I RUN# NTAP E# I T A PES * 
II CELL# IRMIN, IRKAX*ITMIN# ITM^X) 

CALL NEWPAG 

C GO LET USER MAKE A CHOICE AGAIN 

GO TO 33 
C A LA CASA 

50 CALL PIN I TT ( 0# 780 ) 

STOP 

END f ‘ 
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GET START C 
STM 14# 12# 12( 13) 

LR 12,15 
USING GET# 12 
ST 1 3, S AV £*4 
Lft 15,13 
LA 1 3# SAVE 
ST 13,8(,15> 

LM 2 , 6, 0 ( 1 ) 

OC 0(4, 6), 0(6) 

BC 7# ENO 
SW BC 0# LAT ER 
MV 1 SW*l , X • FC * 
OPEN IN 
GET IN 

LM 14,15,20(1) 
SROL 1 4# 4 
SRL 1 5# 4 
ST 1 5# T I ME 
LATER GET IN 
XC 0(4, 6), 0(6) 

LA 7# 32 
LA 6# 4 
LA 8 # 1 
BAL 9,FIXED 
LR 2#3 
LA 3 # 4 
LA 6# 4 
BAL 9# FLOAT 
LM 14, 15, cm 
SRQL 14,4 
SRL 15,4 
S 1 5# T I ME 
ST 15,01,4) 

LA 1,41,1 ) 

LR 2,5 
LA 8,104 
LA 6,28 
BAL 9, FIXED 
RETURN L 13, SAVE ♦4 
LM 14,12,12(13) 

BCR 15,14 
END CLOSE IN 
MV I 3( 6 ), X *FF * 

MV I SW + 1,X*0C* 

BC 15, RETURN 
F IXED LM 14,15,01 1 ) 
SLDL 14,0(6) 

ST 1 4, 0 ( , 2 ) 

LA 2,41,2) 

LA 1 , 4 ( , i ) 

BC T 8, I NC 1 
BCR 15,9 
I NC 1 LA 6, 4 ( , 6 ) 


p 

R 

o 

L 

Q 

G 

U 

E 

GET PARAMETERS 

USER SEDING CLOSE SIGNAL? 

YES, HANOLE IT 
ENTER SWITCH 
SET SWITCH ON 
OPEN FILE FOR INPUT 
GET RUN HEARER RECORD 
PICK UP CLOCK VALUE 
ISOLATE FRACTION BITS 
IN REG 15 

AND SAVE FOR LATER USE 
GET A FRAME 

SIGNAL RECQRO AVAILABLE 
SET SHIFT MAXIMUM AMOUNT 
SHIFT LEFT 4 8ITS 
ONE ITEM 

GO TO FIXED CONVERSION 

SET OUTPUT POINTER 

FOUR ITEMS 

SHIFT LEFT 4 BITS 

GO TO FLOATING POINT CONVERSION 

GET FRAME CLOCK VALUE 

ISOLATE FRACTION BITS 

IN REG 15 

SUBTACT INITIAL FRACTION BITS 
STORE TICKS SINCE RUN STARTED 
BUMP INPUT POINTER 
SET OUTPUT POINTER 
104 ITEMS 

SHIFT LEFT 28 BITS 
GO TO FIXEO CONVERSION 
RETURN 
TO 

USER 

CLOSE FILE 
SET SIGNAL ON 
CLEAR SWITCH 
RETURN 

PICK UP DOUBLE WORD CONTAINING ITEM 

ISOLATE IT 

PUT IT OOWN 

BUMP OUTPUT 

BUMP INPUT 

CONTINUE FOR NUMBER OF ITEMS 
RETURN TO USER 
INCREASE SHIFT BY 4 



_ okGB & 
**$£ <*>*** 


BAL 11 , MORE 
BC 1 5, F I XE D 
FLOAT LM 14,15,0(1) 
SLOL 14,0(6) 

SROL 14,23 
N 14, LAST? 

LR 10,14 
LA 10, 3 ( , 10) 

SRL 10,2 
LR 11, 1C 
LA 11, 64 ( , 1 1 ) 

SLL 10,2 
SR 10,14 
SRL 15,0(10) 

LR 14,11 
SLOL 14,24 
LTR 1 5# 1 5 
BC 11, ♦♦B 
A 14, ONE 
ST 14, 0 ( , 2 ) 

LA 2,4 ( , 2 ) 

LA 1 , 4 ( , 1 ) 

BCT 8, I NC2 
BCR 1 5# 9 
1 NC 2 LA 6, 4 ( , 6 ) 

BAL 11, MORE 
BC 15, FLOAT 
MORE CR 6,7 
BCR 7,11 
LA 1 , 4 ( , 1 ) 

LA 6,0 
BCR 15,11 
SAVE DS 18F 
TIME OS F 

LAST7 DC X* COGOOC7F * 
ONE DC X»00000001* 

IN OCB DSORG*PS, MACR 
END GET 


a? 


AND SET UP FOR MORE 
CARRY OIN 
PICK UP ITEM 
ISOLATE IT 
ISOLATE FRAMTION 
ISOLATE EXPONENT 
SAVE THIS 
BUP TO NEXT 
MULTIPLE OF FOUR 
ANO SAVE THIS 
EXCESS 64 NOTATION 
MULTIPLE OF 4 
GET DIFFERENCE 
AS SHIFT AMOUNT 
PUT THIS ASIDE 
JOIN WITH FRACTION 
ROUND UP? 

NO 

YES, ADD ONE 

AND PUT THE NUMBER DOWN 

BUMP OUTPUT POINTER 

BUMP INPUT POINTER 

CONTINUE FOR ANOTHER ITEM 

OTHERWISE, RETURN TO USER 

INCREASE SHIFT AMOUNT 

CHECK FOR TOO MUCH 

AND CONTINUE 

ARE WE AT 32 BITS? 

NO, CONTINUE 
YES, SKIP ONE FULLWORO 
AND SHIFT 0 BITS 
RETURN TO USER 


F»GL, DDNAME- IN,EOOAD-£ND 





GETS 
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GETS 

START 

0 

P OF POOR *««*“*' 


STM 

14,12,12(13) 

R 


LR 

12,15 

0 


US ING 

GETS, 12 

L 


ST 

13, SA VE*4 

0 


LR 

15,13 

G 


LA 

13, SAVE 

U 


ST 

13, 8( ,15) 

E 


LM 

2, 10, C< 1) 

PICK UP THE PARAMETERS 


OC 

3<4,1C),0( 1C ) 

USER SENDS CLOSE SIGNAL? 


BC 

7, EOD AD 

YES, HANDLE IT 

SW 

BC 

0, L AT ER 

NO, ENTER SWITCH 


MV I 

SW+1, X* FO' 

SET SWITCH ON 


OPEN 

( 01SK,I NPUT) 

OPEN FILE TO BE READ 

LATER 

GET 

0 ISK 

GET A PECORS 


xc 

0(4, 10), 01 10) 

SIGNAL RECORD OBTAINED 


MVC 

0(4,2 ),0(1 ) 

GRID NUMBER 


MVC 

0(4, 3), 4(1) 

FRAME NUMBER 


MVC 

3(4,4 ), 8( 1 ) 

X 


MVC 

0(4,5 ) , 1 2 ( 1 ) 

r 


MVC 

0(4,6 ), 16 ( 1 ) 

R 


MVC 

3(4,71,20(1) 

T 


MVC 

3(4,8 ), 24 ( 1 ) 

CLOCK TICKS 


LA 

1, 28( ,1 ) 

POINT AT ARRAY START 


LA 

10,104 

104 ENTRIES 

AMP 

MVC 

3(1, 9), OC ) 

ONE BYTE EZCH 


XC 

0(3, 9), 0(9) 

CLEZR OUT JUNK 


LA 

l,l(,i) 

BUMP 


LA 

9, 4 ( , 9 ) 

8UMP 


BCT 

10, AMP 

CONTINUE 

RETURN 

L 

13, SA VE+4 

RETURN 


LM 

14,12,12(13) 

TO 


BCR 

15,14 

USER 

EOOAD 

01 

3(10) , X ' FF * 

SET NO RECORD SIGNAL ON 


MV I 

SW*1, x*CO* 

CLEAR SWITCH 


CLOSE 

DISK 

CLOSE FILE 


CLC 

DISK«-40(8),»C , SCRTIN • 

ALTERNATE 


BC 

8, PUT INS 

DDNAMES 


MVC 

D ISK440(8),«C*SCRTIN • 

•SORTIN' 


BC 

15, RETURN 

AND 

PUTINS 

MVC 

D ISK *40 ( 8 ) , *C ' I NS • 

INS 


BC 

15, RETURN 

IN PAIRS 

SAVE 

DS 

18F 


0 ISK 

OC B 

DSORG-PS»MACRF«GL,DDNAME 

« INS, EOOAD«£ODAD,BLK$ I ZE* 12936, 



LRECL •132,R£CFM«FB 



END 

GETS 




SUBROUTINE SORTIT 




SQRTIT START 0 P ' 

STM 14, 12, 12( 13) R 

LR 1 2# 1 5 0 

USING SORT IT# 12 L 

ST 13,SAVE«-4 0 

LR 1 5# 13 G 

LA 1 3, SAVE U 

ST 1 3# 8( # 1 5 ) E 

LINK E P« IERRCOOO,PARAM=I PARM+2 )# VL«1 LINK TO SORTER 

L 13,SAVE*4 RETURN 

LM 14,12,12(13) TO 

BCR 15,14 USER 

DS CF 

PARM DC X'OC’ 

DC AL3 I SORT-SORT A) 

SORTA DC A ( SORT) 

DC A (R ECORD-I ) 

DC A(RECORD) 

DC A( END-1) 

DC A ( C ) 

DC A ( 0 ) 

DC X * F FOCC ID 7 ' 

SORT DC C* SORT F I EL DS * ( 1 , 4 , U I , A ,5, 4, B I , A ), S I ZE » 1 15 100 
RECORD DC C* RECORD TYPE »F, LENGTH® ( 132 ) • 

END ECU ♦ 

SAVE OS 18F 
END SORTIT 
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PUT 

START 

0 


STM 

14* 12 #12(13) 


LR 

12*15 


US ING 

PUT* 1 2 


ST 

1 3* SA VE+4 


LR 

15*13 


LA 

13* SA VE 


ST 

1 3 * 8 ( *15) 


LM 

2* 10*0(1) 

SW 

8C 

0* L AT ER 


MV I 

SW+1* X ' F 0 ' 


OPEN 

( OUT* OUTPUT ) 

LATER 

L 

10*0( *10) 


LTR 

10*10 


BC 

7* C LO S£ 


LA 

1 * BUF FER 


MVC 

0(4*1 )*0(£> 


MVC 

4(4*1 )* C ( 3 ) 


MVC 

8(4*1 ) * 0 ( 4 ) 


MVC 

12(4* 1 ) * 0 ( 5 ) 


MVC 

16(4* 1)*C(6) 


MVC 

20(4* 1 ) *0 ( 7 ) 


MVC 

24(4* 1)*0( 8) 


LA 

l * 28 ( *1 ) 


LA 

10*104 

L CUP 

MVC 

0( 1*1 )*3( 9) 


LA 

1*1(*1) 


LA 

9* 4 ( * 9) 


BC T 

10* L OOP 


PUT 

OUT* BUFFER 

RETURN 

L 

1 3* SA VE +4 


LM 

14,12*12(13) 


BCR 

15*14 

CLOSE 

CLOSE 

OUT 


MV I 

SW*1* x*oo* 


BC 

1 5* RE TURN 

SAVE 

DS 

18F 

BUFFER 

DS 

33 F 

OUT 

DC B 

DSORG*PS*DDNAME 


BLKSI Z£ *12936 
PUT 


P 

ft 

0 

L 

□ 

G 

U 

E 

PICK UP PARAMETERS 

ENTER SWITCH 

SET SWITCH ON 

OPEN FILE TO BE WRITTEN 

USER SENDS 

CLOSE S IGNAL? 

YE S* HANDLE IT 
POINT TO RECORO AREA 
GRID NUMBER 
FRAME NUMBER 
X 
Y 
R 
T 

r i nr k t i r k s 

POINT AT ARRAY OUTPUT AREA 

104 ELEMENTS 

TAKE ONE 

BUMP 

BUMP 

CONTINUE 

WRITE OUT RECORD 

RETURN 

TO 

USER 

CLOSE FILE 
CLEAR SWITCH 
RETURN 


SORTIN*MACRF*PM*RECFM*FB*LR£CL«132# 


END 



SOs *<Wriiffi 


3s 


digits 


Q 
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SUBROUTINE 0 IG ITS C JVAL, L INE, N ) 

SUBROUTINE TO PRODUCE TEKTRONIX PRINTABLE OIGITS FROM A NUMBER 
JVAL IS INTEGER VALUE WITH NO MORE THAN 6 OIGITS. POSITIVE. 
LINE IS THE ARRAY TO RECEIVE RESULTS 
N IS THE NUMBER OF DIGITS DESIRED 
DIMENSION L INE (6) 

C FORM FIRST DIV ISOR 

I NDEX-IO** ( N-l ) 

C LOOP ONCE FOR EACH DIGIT 

DO 1 I «1 > M 

C GET QUOTIENT 

IF- JVAL/ INDEX 

C NUMBER REDUCED TO REMAINDER 

JVAL-JVAL-IF*I NDEX 

C WE HAVE PRINTABLE TEKTRONIX INTEGER 

LINE! I ) * IF *48 

C REDUCE DIVISOR FOR NEXT TIME 

1 INDEX- INDEX/LO 

C ALL DONE 

RETURN 
END 


ORIGINAL PAGE IS . 
OF POOR QUALflt 
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r>r>f>r, or>r>r> 
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***♦ *************** *************************** ****** ** *************** 

* 

SUBROUTINE MENU * 

* 


THIS SUBROUTINE GIVES USER A CHOISt OF "SCAN'* DISPLAY * 
(RADIUS X ANGLE X EITHER VELOCITY OR AMPLITUDE) OR * 
"RAY" DISPLAY (RADIUS X VELOCITY X AMPLITUDE). * 


C 


SUBROUTINE MEN U ( ISTART, Ti, ISTOP, T2, ISCAN, ISCANi, IDAY, 

1 IRUN, NT APE , ITAPES, (CELL* I RM IN* IRMAX, 

2 ITMIN* ITMAX ) 

DIMENSION ISTART (60)* Tl(60)* IST0P(60)* T2<60), ITAPES(2) 

MTAP E«NT APE 

1 CALL MOVABS ( 80* 690 ) 

C DISPLAY CHO IS ES 

CALL AOUTS T( 7* ’OPTIONS* ) 

CALL MOVABS ( 100* 660 ) 

CALL AOUTS T ( 45* *1 ... THREE DIMENSIONAL ONE FULL SCAN DISPLAY’) 
CALL M 0 V A 8 S ( 100* 630 ) 

CALL A 3UTS T( 39* *2 ... THREE DIMENSIONAL ONE RAY DISPLAY’ ) 

CALL MOVA BS ( ICC, 6CC ) 

CALL AOUTST( 15* '9 ... TERMINATE’ ) 

C ENTER OPTION 

90 CALL A I NST ( 1* IGPT » 

CALL C ORE ( I3PT* 1 ) 

READ ( 99* 901 ) IOPTN 

90 1 FORMAT ( 11 ) 

IF ( IOPTN .EQ. 9 ) GO TO 999 
IF ( IOPTN .NE. i ) GO TO 2CC 
C USER CHOSE "SCAN" DISPLAY OPTION 

CALL AOUTST ( 36* ’ E NTE R STARTING SCAN NUMBER (2 OIGITS)’) 

CALL A INST( 2* I SC N) 

CALL C ORE ( I SCN *2 ) 

READ (99,902) I SC N 
902 FORMAT (12) 

DO 120 NFRM*ISCN» ISCAN 

IF (MTAPE.E0.2 .AND. NFRM. GE . I SC AN1 ) NT AP E *3 

CALL D I SPL Y( ISTART* Tl, I S TOP* T2» I SC AN* ISCAN1* IDAY* IRUN* 

1 NTAPE * ITAPES* ICELL* NFRM* ITMIN* ITMAX ) 


120 CONTINUE 


GO TO 1 

C OPTION MUST BE 1* 2 OR 9 

200 IF ( IOPTN .NE. 2 ) GO TO 9C 
C USER CHOSE M R AY” DISPLAY OPTION 

520 CALL ERASE 

CALL MOVABS ( IOC* 750 ) 

CALL AOUTS T ( 12* ’ENTER SCAN #• ) 

CALL A I NST( 2* NPRM ) 

CALL CORE ( NFRM * 2 ) 

READ ( 99, 902 ) NFRM 

IF ( NFRM .GT. NSC AN ) GO TO 52C 

CALL ERASE 


0 ? * 
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CALL RAY( ISTART, Tl, ISTOP* TZ, I SCAN* I SC ANl# IOAY* IRJN* 

1 NT AP t * ITAPtS. ICELL* NF KM# ITMIN* UMAX ) 

GO TO i 

999 CALL F I N I T T ( C# 0 ) 

STOP 

END 
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**»*******#**«*«*** ************************************************** 

* 

SUBROUTINE DISPLAY * 

* 

THIS SUBROUTINE OISPLAYS ONE "SCAN” OF DATA IN THREE * 

DIMENSIONAL form * 

* 


SUBROUTINE DISPLYI NFRMi# ANGI# NFRM2# A N G 2 # NS CAN# ISCAN1# IDAY# 

IRUN# NTAPE# ITAPES# ICELL# ISC AN# ITM IN# ITMAX I 
NFRMKtQ)# ANG1I6G)# NF RM2 ( 60 ) # ANG2I60)# ITAPES(2) 
XB(ICOO)# XWtlGCC) 

XV ( LOGO > # NF(IOO) 

IaMPUOa)# X Y ( 4 ) # IN I 2 ) 

XX(IOCO)# XA ( 1 OOC ) # N X ( 1 00 ) » TA3L E ( 100 ) 

ZZl 100 ) 

Z (ICO# ICC I# NGRIDI4) 

XA1( 10D0)#XVi< 1000)#XW1( 1000) 

HARP / IC AL# NSCOPE 

# AN GA V E ( 50 ) 


9o2 


963 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
COMMON / 

DIMENSION ANGI 1000) 

DIMENSION I I ( 5 C ) 

DOUBLEPRtC IS ION LFRM1# AANG1# LFRM2# AANG2# LTAPE1# LTAPE2 
DOUBLEPRECISION ANGAV ( 5C ) # TABL E2 ( 50 ) 

DOUBLEPRcC ISION ZAVE2# ZH IGH2 
NT A PE 3 » 3 
ISAME « 0 

ANGLES GIVEN ON NASA TAPE ARE ONE HALF OF ACTUAL ANGLE S 
A ANGI « ANGl(ISCAN) * 2.0 
AANG2 • ANG2IISCAN) * 2.C 

CONVERT NUMERICAL FORM TO CHARACTER FORM FOR DISPLAY 
WRITE ( 2# 962 ) NF RK 1 < I SCAN ) # A ANGI # NFRM2 ( I SC AN ) # 

L AANG2# IDAY# IRUN, ITAPES# ISCAN# ICf 

FORMAT ( 16# Ft. 3# 16# Fb.3# 13# 12# 215# 12# 13 
REWIND 2 

READ ( 2# 9t3 ) L F R M 1 # AANG1# LFRM2# 

L LTAPE2# LSCAN# LC ELL 

FORMAT ( 4A&# A3# A2# 2a5# A2# A3 ) 

REWIND 2 
CALL ERASE 

DISPLAY INFORMATIONS FOR THIS RUN 


) 


AANG2# LDAY# LRUN# LTAPEi# 


CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
IF 
IF 


CHRSIZI2 ) 


MOVABSI 
AOUTSTI 
AOUTSTI 
AOUTST ( 
AOUTSTI 


1G0# 750 ) 

32# 'NASA PROJECT DUST DEVIL 


DAY 


» 


RUN # • ) 


IF 
IF 
IF 
CALL 
CALL 


3# LDAY 
10 # * 

2# LRUN ) 

INTAPE.EQ .3 ) CALL A0UT$T(16#' TAPE 2 GF 2-') 
INTAPE.EQ. 2 ) CALL AOUTST (16#* TaPE 1 OF 2-» 
I NTAPE .EQ. 1 ) CALL AOUTSTI 13# • TAPE # 

INTAPE • EQ. 3) CALL AOUTST !5#LTAPE2) 

INTAPS. LE. 2) CALL AOUTST !5#LTAPE1> 

C HRS I Z I 3 ) 

MOV A3 S I 30# 725 ) 


) 

1 -' 



CALL AGUTS T ( 7* ’SCAN » • ) 

CALL AOUTSTI 2* LSCAN ) 

CALL AOUTSTI 13* • FRAME #'• ) 

CALL AGUTS T ( 6 * LFRMi ) 

CALL AGUTS T ( 9, • ANGLE • ) 

CALL AOUTSTI t, A A NG 1 ) 

CALL AuUTSTI 18* ' TO FRAME U • ) 

CALL ACUTSTI 6* LFRM2 ) 

CALL AOUTSTI 9* '• ANGLE • ) 

CALL AOUTSTI 6, A A NG2 ) 

CALL AOUTSTI 25* ' CRITICAL AMP AT CELL • ) 

CALL AOUTST I 3* LCELL ) 

CALL MOVABSI 8C*b9C ) 

DISPLAY OPTIONS FOR UNTRANSLATED DATA 
CALL AOUTST I 7* 'OPTIONS* J 
CALL MGVABS I ICO* 660 ) 

IF I IOAY .GE. 238 ) GO TO 605 
CALL AOUTSTI 3* *11 ... I* ) 

CALL CHRSIZI4I 

CALL AOUTSTI 3* 'MAX' ) 

CALL CHRSIZI3) 

CALL AOUTSTI 15* '* R* A DISPLAY' ) 

CALL MOVABSI 100* 63C ) 

CALL AOUTSTI 37* '12 ... I* R* A DISPLAY AT GIVEN V* 
CALL MOVABS I 100* bOO ) 

CALL AOUTSTI 39* '21 ... V* R, A OISPLAY (VELOCITY AT 
CALL CHRSIZU) 

CALL AOUTSTI 3* 'MAX' ) 

CALL CHRSIZI3) 

CALL AOUTSTI 1* 1 

CALL MOVABS I ISO* 570 > 

CALL AOUTSTI 8* '22 ... V' > 

CALL CHRSIZI4) 

CALL AOUTSTI 3* 'MAX' ) 

CALL CHRSIZI3) 

CALL AGUTSn 13* ' R* A DISPLAY* » 

CALL MOVABS I ICO* 540 ) 

CALL AOUTSTI 12* *31 ... V * I' ) 

CALL CHRSIZI4) 

CALL AOUTSTI 3* 'MAX' ) 

CALL CHRSIZI3) 


CALL 

AOUTSTI 

1 5* ' 

* R * A 

DISPLAY' ) 

CALL 

MOVABS 

1 100 

*. 510 

) 

CALL 

AOUTSTI 

2 C* • 55 ... 

"RAY" DISPLAY' ) 

CALL 

MOVABS 

I IOC 

* 4 8C 

) 

CALL 

AOUTSTI 

22* ' 

99 ... 

GO TO NEXT DATA 

CALL 

C HRS I Z ( 4 > 



CALL 

MOVABSI 

100* 

450 ) 


CALL 

AOUTST 

1 15* 

•MAX 

... MAXIMUM' ) 

CALL 

MOVABSI 

100* 

425 1 


CALL 

AOUTST 

1 15* 

•MIN 

... MINIMUM' ) 

CALL 

MOVABSI 

ICO* 

4C0 ) 


CALL 

AOUTST 

( 15* 

' I .. 

. INTENSITY* ) 

CALL 

MOVABSI 

100* 

375 ) 
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) 


I* 
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505 

C 


C 

6C 5 


CALL 

A OUTS T I 14# 

•V . 

.. VELOCITY* 

CALL 

M0VA8SI IOC# 

350 

) 

CALL 

AOUTS T I 12# 

•k . 

.. RADIUS* ) 

CALL 

MCVABSI 100# 

32 5 

) 

CALL 

AOUTS T I 11# 

* A . 

.. ANGLE' 1 

CALL 

CHRSIZI3 I 



CALL 

M0VA8SI 100# 

28 0 

) 

CALL 

AOUTSTI 13# * 

ENTER 

OPTION* I 


ENTER UPTIQN FOR UNTRANSLATED DATA 
CALL A I NS T ( 2# IQPTN ) 

CALL CORE ( IGPTN# 2 > 

READ ( 99, 90 2 ) IQPTN 

IF ( IQPTN .EQ. 55 > GO TO 710 

IOP*IOPTN 

GO TO 5Co 

DISPLAY OPTIONS FOR TRANSLATED DATA 
CALL AOUTSTI 9, 'ill ... I* ) 

CALL CHRSIZ14) 

CALL AOUTSTI 3# ’MAX' I 
CALL CHRSIZI3I 

CALL AOUTSTI 4 2 , •# R# A OISPLAY (PART i OF 2# CELL 5 2 
CALL MOVABSI 100* 630 ) 

CALL AOUTSTI 43/ *112 ... PART 2 OF OPTION 111# CELL 53 ■ 
CALL MCVA6S I IDG# 600 ) 

CALL AOUTSTI 44#*121 ... I# R# A DISPLAY AT GIVEN V 
CALL MOVABS I 100# 5 7C ) 

CALL AOUTSTI 40#*2U ... V# R# A DISPLAY (VELOCITY AT 
CALL CHRSIZI4) 

CALL AOUTSTI 31# 'MAX. PART 1 OF 2# CELL 52 - 1 • ) 

CALL CHRSIZI3I 

CALL AOUTSTI 1 # • > • I 

CALL MOVABSI 100# 540 ) 

CALL AOUTSTI 8# *212 ... • ) 

CALL CHRSIZI4) 

CALL AOUTSTI 35# 'PART 2 OF OPTION 211# CELL 53 - 104* ) 
CALL CHRSIZI3) 

CALL MOVABSI 100# 510 I 
CALL AOUTSTI 9#*22i ... V* ) 

CALL CHRSIZI4) 

CALL AOUTSTI 3# ‘MAX' ) 

CALL CHRSIZI3) 

CALL AOUTSTI 13# * R# A DISPLAY* ) 

CALL CHRSIZI4) 

CALL AOUTS T 1 2 9 # • (PART 1 OF 2# CELL 52 - 1 I* I 
CALL CHRSIZI3) 

CALL MOVABSI ICC# 480 ) 

CALL AOUTSTI 8# * 222 ... * ) 

CALL C HRS I Z I 4 ) 

CALL AOUTS T( 3 5 # ’PART 2 OF OPTION 221# CELL 53 - 1C4* ) 
CALL CHRSIZI3) 

CALL MOVABSI 100# 450 ) 

CALL AOUTSTI 1 3# • 31 1 ... V * I ' ) 

CALL C H R S I Z I 4 ) 

CALL AOUTSTI 3# 'MAX' ) 


- II* I 
104* ) 
CELL* I 
I* 1 
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CALL CHRSIZI3) 

CALL AOUT$T( 15# *# R# A DISPLAY* ) 

CALL CHRSIZ(A) 

CALL AOUTSTI 29# • (PART 1 OF 2# CELL 52 - 1 ) ' ) 

CALL CHRSIZ(i) 

CALL MOV A BS ( 100# 420 ) 

CALL AOUTSK 8# *312 ... * ) 

CALL C HRS I Z ( 4 ) 

CALL AOUTST ( 35 # 'PART 2 OF OPTION 311# CELL 53 - 104* I 
CALL MOVA 6S ( 100# 3bO ) 

CALL CHRS I Z ( 3 ) 

CALL AOuTS T ( 28# ' 555 OR 55 ... "RAY” OISPLAY* ) 

CALL MCVABSI 100# 350 ) 

CALL AOUTS T 1 2 9 # * 999 OR 99 ... GO TO NEXT DATA* ) 

CALL MOVA BS ( 100# 300 ) 

CALL AOUT$T( 13# 'ENTER OPTION* ) 

C ENTER OPTION FOR TRANSLATED DATA 

CALL A I NS T ( 3# IOPTN ) 

CALL CORE ( IOPTN# 2 > 

READ < 99 , 903 ) IOPTN 
I 0P» I0PTN/1C 


IF 

( 

IOPTN .EQ. 

55 ) 

GO 

TO 

710 

IF 

( 

IOPTN .EQ . 

550 

) GO 

TO 

710 

IF 

( 

IOPTN .EQ. 

555 

) GO 

TO 

710 


50b I F ( I0P.NE.22 ) GO TO 507 
CALL MOVA BS ( iGO# 260 ) 

CALL AOUTST ( 43# 'ENTER AMP FOR MAX VELOC ITY . .. ( 3 DIGITS)' ) 

C ENTER AMPLITUDE THRESHOLD 

CALL A I N S T ( 3# IVAMP ) 

CALL CORE ( IVAMP# 2 ) 

READ ( 99# 903 ) IVAMP 
5*- 7 CONTINUE 

C IF OPTION X IS 99 OK 999 RETURN TO CALLING PROGRAM 

IF ( IOPTN .cQ. 99 ) RcTU.RN 
IF ( IGPTN.E0.999 ) RETURN 
IF ( IOPTN. E 0.9 9C ) RETURN 
CALL ERAS £ 

C SELECT TAPE IF MORE THAT ONE TAPE USED 

IF < MFRM2 .Nt. N FR Mi ( I S CAN ) ) GO TO 512 
IF ( NTAPE .EQ. 3 .AND. NTAPE3 .EQ. 0 ) GO TO 512 
XX<1> * XX(NN) 

ANG(l) * ANG(NN) 

XV(i) » XV(NN) 

X A ( 1 ) - XA(NN) 

X B ( 1 ) « XB(NN) 

XW(1) * XW(NN) 

XA1 ( 1 ) * XAi < NN ) 

X B ( 1 ) * XB ( NN ) 

XVH 1 ) «XVi ( NN ) 

XW1 ( 1 ) *XWl ( NN) 

GO TO 5.6 

C READ DATA 

C IF NT APE * 1 OR 2 USE FIRST TAPE 

512 IF ( NTAPE -L£. 2 ) CALL GET (IFRM# XY# ITM# IAMP# IRET) 
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C IF NPATE *3 USE SECOND TAPE 

IF I NTAPE .EQ. 3 I CALL GET21IFP.M, XT, ITM* IAMP* IRET) 

IF I NTAPE .EQ. 3 ) NTAPE3 « 1 
IF ( IRET .NE. 0 ) IRET » 0 
IF I NFRMit ISCAN) .NE. IFRM ) GO TO 512 
C STORE DATA FOR LATER USE 

WRITE ( 3* <964 ) IFRM* XY14)* IAMP 
964 FORMAT ( 16* F5.2* 1C4I3 ) 

XXII) * X Y I 3 ) 

ANG(I) « XYIO * 2.C 
IF! IDAY.Gh.238 ) GO TO 610 
MAMP « I AMP ( 1 ) 

XVII) * 1 

C FIND MAXIMUM AMPLITUDE AND VELOCITY 

DO 514 I * 2* IOC 

IF I I AMP ( I ) .GT. MAMP ) XV(i) - I 

IF I I AMP ( I ) .GT. MAMP ) MAMP « I AMP I I ) 

514 CONTINUE 

X A ( 1 ) ■ MAMP 
X 6 1 1 ) • IAMPI I CELL ) 

XWC) * L 

C FIND MAXIMUM VELOCITY USING AMPLITUDE THRESHOLD FOR 

C UNTRANSLATED DATA 

DO 515 I » 2* 1GG 

IF I I4MPII) .GT. IVAMP 1 XW(l) • I 

515 CONTINUE 
GO TO 516 

610 X V i 1 1 ) ■ A . j 
X V I 1 ) * I . G 
MAMP-IAMPI52) 

MAM PI * IAMP I 53) 

J B 5 1 

C FIND MAXIMUM AMPLITUDE AND VELOSITY FOR TRANSLATED DATA 

DO 612 1*54* iOC 
IFC IAMPI I J.GT.MAMPl ) XV1U)*I 
IF! IAMPI I ).GT. MAMP1 ) MaMPI* IAMPI I ) 

IF I IAMPI J ).GT. MAMP > XVI l ) » I 

IF I IAMPI J ) .GT .MAMP ) MAMP-IAMPIJ) 

J-J-i 

oi 2 CONTINUE 
X A I 1 ) * MA MP 
XAl li )*MAMP: 

X 8 1 1 ) * I A MP I ICE LL ) 

XWID-1.0 

xwiiD-i.o 

J»51 

C FIND MAXIMUM VELOCITY FOR TRANSLATED DATA 

DO 614 I = 5 4* 1 0 C 
IF! IAMPI I ).GT. IVAMP ) XWl(i)*I 
IF! IAMPI J).GT. IVAMP) XWli)*I-52 
J* J-i 

6i 4 CONTINUE 

516 MFRM2 - NFRM2I ISCAN) 

NN = NFR M2 1 ISCAN) - NFRMUISCAN) ♦ 1 
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DO 523 I - NN 
C READ DATA 

IF ( NT A P £ .LiE. 2 ) CALL GET (IFfifW XT# I TM# IAMP# IRET) 

IF ( NTAPE . E Q . 3 > CALL G2T2IIFRM# XT# I TM# IAMP# IRET) 

C STORE DATA FOR LATEk USE 

WRITE ( 3# 964 ) IFRM# XTI4)# IAMP 

X X ( I ) * XT ( 3 ) 

ANG ( I ) = X T I 4 ) * 2.C 

IF ( IDAT.GE.238) GO TO 62C 

C FOR UNTRANSLATED DATA 

MAMP • IAMP(I) 

XVIII » 1 

FIND VELOCITT AND AMPLITUDE FOR UNTRANSLATED DATA 
DO 522 J ■ 2# 100 

IF ( IAMPIJ) .GT. MAMP ) XV(I> « J 

IF ( IAMPIJ) .GT. MAMP ) MAMP * IAMPIJ) 

522 CONTINUE 

X A ( I ) * MAMP 
X B ( I ) * IAMP(ICELL) 

X W ( I > > x 

FIND MAXIMUM VELOCITT USING AMPLITUDE THRESHOLD FOR 
UNTRANSLATED DATA 
DO 524 J - 2# ICC 

IF I IAMPU) .GT. IVAMP ) XWII) * J 
524 CONTINUE 
GO TO 5<.8 

C FOR TRANSLATED DATA 

62 G MAMP- IAMP! 52 ) 

MAMP1-IAMP (53) 

XVI I) *1.0 
XVI ( !)■:.£> 

K-51 

C FIND MAXIMUM POSITIVE VELOCITT 

DO 622 J»54#10C 

I F ( IAMPIJ). GT.MAMP1) XVI I I )* J-52 
IF! IAMPIJ ).GT. MAMP 1) MAMP1 * I AMP I J ) 

C FIND MAXIMUM NEGATIVE VELOCITT 

IF! IAMPIK). GT. MAMP) XV(U*J-52 
IF I IAMPIK) .GT .MAMP) MAMP « I AMPI K ) 

K-K-l 

622 CONTINUE 

XAI I ) * MAMP 

XA1 I I ) * MAM P 1 

XBI I ) • I A M P I ICE LL ) 

XWI I 1*1.0 
XW1UW.9 
K* 5 1 

C FIND MAXIMUM POSITIVE VELOCITT USING AMPLITUDE THRESHOLD 

DO 62 4 J * 5 4# 1 0 G 

IF! IAMPIJ ).GT. IVAMP) XWl(I)«J-52 

C FIND MAXIMUM NEGATIVE VELOCITT USING AMPLITUDE THRESHOLD 

IF! IAMPIK). GT. IVAMP ) XW(I)*J-52 
K*K-1 

624 CONTINUE 
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52o 

C 

53l 

534 

536 

540 

542 

546 

5*»8 


CONTINUE 
REWIND 3 

FIND GRID SUE FOR BASE OF 3-D DISPLAY 
KNT = 2 
N X ( i ) * _ 

NF ( i ) = NFRMK ISCAN ) 


I » 1 

IF ( X X ( I > .GT. X X ( 2 > ) GO TO 534 

I * I ♦ • 

IF ( H-i .GE. NN ) GO TO 536 

IF ( XX(I> .LE. XX(I*1) ) GO TO 532 

NX! KNT ) * I 

NF(KNT) « NFRMi(ISCAN) + I 
KNT » KNT ♦ 1 
I » I + i 

IF ( 141 . GE . NN ) GO TO 536 

if ( xxm .gt. xx<i+n > go to 534 

NX(KNT) • I 

NFIKNT) ■ NF RM i ( ISCAN) ♦ I 
KNT « KNT ♦ i 
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GO TO 532 

NX(KNT) ■ NN 

NF(KNT) * NF RM2 ( ISCAN) 

SUMA « 0.0 

DO 540 l«l> KNT, > 

S UNA » SUMA ♦ XX ( NX ( I ) ) 
SUMB » 0.0 

DO 542 I « 2# K NT * 2 
SUMB « SUMB «■ XXINXI 1)1 
C » ( KNT+i ) t 2 
ASUHA * SUMA / C 
C « KNT / 2 
ASUM8 * SUMB / C 


IF 

( 

ASUMA .GT. ASUMB 

) 

AMAX 

« ASUMA 

IF 

( 

ASUMA .GT. ASUMB 

) 

AMIN 

« ASUMB 

IF 

( 

ASUMA .LT. ASUMB 

I 

AMIN 

« ASUMA 

I F 

( 

ASUMA .LT. ASUMB 

) 

AMAX 

« ASUMB 


I A VE * NN / KNT 
TABLE ( 1 ) * AMIN 
T ABLE I IA VE ) • AMAX 
MM * IAV6 - 1 
DO 546 I ■ 2# MM 

TABLE! I » • TABLE(I-i) ♦ ( AMAX - AMIN > / FLOAT ( IAVE ) 

KY - IAVE 

KX - KNT • l 

DO 556 I « l# KX 

DO 554 K » 1, KY 

XXMIN « iOCO.O 

NNX * Nxm 

NNNX * NXUol) 

DO 552 J » NNX, NNNX 

XXM * ABS( TABLE(K) - XX(J) ) 

IF ( XXM .GE. XXMIN ) GO TO 552 
XXMIN = XXM 
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jxx « J 
552 CONTINUE 

C TRANSFER DATA FOR 3-D DISPLAY 


IF 

I 

IOPTN .EU . 11 

) 

Z< I*K) 

m 

XAI JXX) 

IF 

I 

IOPTN .EQ. 12 

) 

Z( I*K) 

8 

XBI JXX) 

IF 

I 

IOPTN .EQ. 21 

) 

Zt I*K) 

• 

XVI JXX) 

IF 

I 

IOPTN .EQ. 22 

) 

ZI I * K ) 

8 

XWI JXX) 

IF 

I 

IOPTN .EQ. 21 

) 

Zl I*K) 

m 

XAIJXX) * XVIJXX) 


IF( IOPTN.EQ.L11 ) Zt I*K)»XA( JXX) 

IF( IOPTN.EQ.112) Z I I # K ) »XAi I J XX ) 

IF( I0PTN.EQ.121) Z< I*K)»XBIJXX) 

IFdOPTN.EQ.2U) Z ( I * K ) » XV ( J XX ) 

I FI I0PTN.EQ.212) Z d * K ) » XVI I J XX ) 

IF! IQPTN.EQ.il 1) Z ( I * K ) « XA ( J XX > *XV I J XX ) 

IFI I0PTN.EQ.312) Z 1 I * K ) « XA i ( J XX ) *X VI ( J XX ) 

IF! I0PTN.E0.221) Z I I * K ) ■ XW I JXX ) 

IFI I0PTN.EQ.222> Z 1 I * K ) • XW1 1 J XX ) 

55* CONTINUE 
556 CONTINUE 

C FIND SCALING FACTOP FOR 3-D DISPLAY 

ZLOW » ZI1*1) 

ZHIGH • ZU*l) 

DO 558 I * l» KX 
DO 558 K«i* KY 

IF I Z I I # K ) .GT. ZHIGH ) ZHIGH » ZII»K) 

IF I Z I I > K ) .LT. ZLOW ) ZLOW - ZU#K) 

558 CONTINUE 

Z SUM » 0.0 

C SCALE DATA FOR 3-D DISPLAY 

DO 560 I ■ 1 » KX 
DO 560 K » 1/ KY 
Z SUM « ZSUM ♦ Z I I * K ) 

Z I I > K ) • I Z I I > K ) - ZLOW ) / I ZHIGH - ZLOW I * 200.0 
560 CONTINUE 

ZAVE * ZSUM / FLOAT! KX ♦ KY ) 

IF I ANG1IISCAN) .IE. AMG2I ISCAN) ) GO TO 566 
C IF DATA IS PLACED BACKWARD* REPLACE DATA 

DO 565 K * 1* KY 
DO 562 I « 1* KX 
562 ZZII) * Z I I * K ) 

KK • KX 

DO 56* I - 1» KX 
ZI I*K) » ZZIKK ) 

KK * KK - 1 
56* CONTINUE 

565 CONTINUE 

C SET THE OPTIONS FUR 3-D DISPLAY 

566 NPY « KX 
NPX « KY 
NBORD - 0 
NGRIDU) - 0 
NCR 10(2) • 0 
NGRIDI3) * 0 
NGR I D I * ) » C 
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ICAL * 0 
ZMAX » 30.0 
R = 200. C 
THETA » 80.0 
PHI « 30.0 
Z MI N = 0.0 
NHIDE * 1 
NSC AL E « 2 

IF ( I0PTN .EQ. 21 .OR. IOPTN .EQ. 22 ) ZMAX- 15.0 
CALL CHRSIZI2) 

C LABEL DISPLAY 

CALL MOVABS I ICO* 759 ) 

CALL AOUTSTI 32* * NA SA PROJECT OUST OEVIL OAY • I 

CALL AOUTS T ( 3* LOAY ) 

CALL AOUTS T I 10* ' RUN f ' ) 

CALL AOUTSTI 2* LRUN ) 

IF (NTAPE.EQ.3 ) CALL AOUTST ( 16* • TAPE 2 OF 2-' I 
IF (NTAPE.EQ.2 ) CALL AOUTST 116*' TAPE i OF 2-') 

IF ( NTAPE .EQ. 1 ) CALL AOUTSTI 13* • TAPE « 1-' ) 
IF INTAPE .EQ.3) CALL AOUTST I5*LTAPE2) 

IF INTAPE. LE. 2) CALL AOUTST I5*LTAPE1I 
CALL CHRSIZI3) 

CALL MOVABS I 7C* 740 ) 

CALL AOUTSTI 7* 'SCAN # • > 

CALL AOUTSTI 2* LSC AN I 

CALL AOUTSTI 13* • FRAME t • ) 

CALL AOUTSTI 6* LFRM1 ) 

CALL AOUTSTI 9* • ANGLE • ) 

CALL AOUTSTI 6* AANG1 ) 

CALL AOUTSTI 18* • TO FRAME # • ) 

CALL AOUTSTI 6* LFRM2 ) 

CALL AOUTSTI 9* • ANGLE • > 

CALL AOUTSTI 6* AANG2 ) 

CALL AOUTSTI 25* • CRITICAL AMP AT CELL ' I 

CALL AOUTST I 3* LCELL I 

CALL MOVABS I 70* 71G) 

CALL AOUTSTI 8* 'OPTION #• > 


IF 

I 

IOPTN .EQ. 

11 

> 

CALL 

AOUTSTI 

3* 

• 

11' 

IF 

1 

IOPTN .EQ. 

12 

) 

CALL 

AOUTSTI 

3* 

• 

12' 

IF 

1 

IOPTN .EQ. 

21 

) 

CALL 

AOUTSTI 

3* 

• 

21' 

IF 

1 

IOPTN .EQ. 

22 

) 

CALL 

AOUTSTI 

3* 

• 

22* 

IF 

1 

IOPTN .EQ. 

31 

) 

CALL 

AOUTSTI 

3* 

9 

31* 


IFI IOPTN. EQ. Ill) CALL AOUTSTI 5* ' 11-1' 
IFI IOPTN. EQ. 112) CALL AOUTSTI 5* ' 11-2' 
IFI IOPTN. EQ. 121) CALL AOUTSTI 5* ' 12-1' 
IFI IOPTN. EQ.211) CALL AOUTST 15*' 21-1* 
IFI IOPTN. EQ. 212) CALL AOUTST 15*' 21-2' 
IFI IOPTN. EQ. 221) CALL AOUTST I 5* • 22-1* 
IFI IOPTN. EQ. 222) CALL AOUTST 15*' 22-2' 
IFI IOPTN. EQ. 311) CALL AOUTST I 5* ' 31-1' 
IFI IOPTN. EQ. 312) CALL AOUTST I 5 *• 31-2 ' 
IFI IOPTN. GT. 300) GO TO 567 
IF I IOPTN .EQ. 31 ) GO TO 567 
WRITE I 2* 969 ) ZHIGH* ZAVE 



969 FORMAT ( 2F 6 . 1 ) 

REWIND 2 

C TRANSLATE real number to character form 

READ ( 2, 97 U ) ZHIGH2, ZAVE2 

970 FORMAT I 2X, A6, 2X, A6 ) 

REWIND 2 

CALL MOVABS ( C, 675 ) 

CALL A OUTS T ( 7* ‘MAXIMUM* ) 

CALL CHRSIZI2) 

IF! IQPTN.GT.10C) IOP« I0PTN/10 

C DISPLAY MAXIMUM AND AVERAGE AMPLITUDE JR VELOCITY 

IF ( IOP .LT. 20 ) CALL AOUTSTI 3, * I ' ) 

IF ( IOP .GE. 20 > CmLL AOUTSTI 3, * V ‘ I 

CALL CHRSIZI3) 

CALL AOUTSTI 1, •*• ) 

CALL AOUTSTI 6, ZHIGH2 ) 

CALL MOVABSI C, 65C > 

CALL AOuTS T I 7, ‘AVERAGE* ) 

CALL C HRS 1212 ) 

IF I IOP .LT. 20 ) CALL AOUTSTI 3, ' I * I 

IF I IOP .GE . 20 ) CALL AOUTSTI 3t * V * I 

CALL CHRSIZI3) 

CALL AOUTSTI 1, •*• ) 

CALL AOUTSTI bp ZAVEZ ) 

567 CONTINUE 

IF I I SAME .£0.1 > GO TO 572 
KNTi » KNT-1 
DO 570 I - 1, KNTI 
K • NX! I ) 

L « NX! 1+1 ) 

KL • C 

ANGAVEI i ) - 0. G 

DO 568 J » Kp L 

ANGAVE ( I ) • ANGIJ) + ANGAV £ ( I ) 

KL * KL ♦ 1 

568 CONTINUE 

ANGAVEII) » ANCAVEII) / FLOAT I KL ) 

57 J CONTINUE 
ISAME = 1 

C TRANSLATE NUMERICAL DATA TO CHARACTER FORM 

WRITE I 2, 968 ) I Ip 1*1, 50) 

968 FORMAT I 2CI I Zp 2X ) ) 

REWIND 2 

READ I Zp 966 ) I I 
966 FORMAT I 2&AA ) 

REWIND 2 

WRITE I Zp 967 ) I ANGAVEI I ) p 1*1, KNTi ) 

9o 7 FORMAT I ICF8.1 ) 

REWIND 2 

READ I Zp 972 ) I ANGAV I I ) * 1*1, KNTI ) 

972 FORMAT I 101 AX, A A ) ) 

REWIND 2 

WRITE I 2, 967 ) I TABLE I I ), 1*1 , I AVE ) 

REWIND 2 
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READ ! 2# 973 ) ( TABLE2 < I )# I»1#IAV£ > 

973 FORMAT ( 1 0 ( 3 X, A5 )) 

REWIND 2 

572 IY - 575 

C DISPLAY ANGLES 

CALL CHRSIZI3) 

CALL MOVABSC 15 * 59C ) 

CALL AOUTSK 5# 'ANGLE* ) 

CALL CHRSIZ(4> 

DO 573 I « 1# KNT 1 
CALL MOVA BS ( 5* IY ) 

CALL AOUTSK 4# 11(11 ) 

CALL AOUTS T ( 4, ANGAV< I ) ) 

IF I KNT1 .GE. 40 I IY ■ IY - 14 

IF ( KNT1 .LT. 40 ) IY * IY - 15 

573 CONTINUE 

CALL MOVABS ( 68C# 4C ) 

IF ( ANGKISCAN) .LE. ANG2 ( I SC AN ) ) CALL AOUTST< 

IF ( ANGKISCAN! .GT. ANG2 ( I SC AN ) ) CALL AOUTSK 

CALL MOVABS ( 180# 560 ) 

IF ( ANGKISCAN) .LE. ANG2 ( I SCAN ) ) CALL AOUTSK 
IF ( ANGKISCAN) .GT. ANG2USCAN) ) CALL AOUTSK 
CALL CHRS I Z( 3 ) 

CALL MOVA BS ( 960# 615 ) 

CALL AOUTS T ( 1# 'R' ) 

IY » 60C 

CALL CHRSIZ44) 

C DISPLAY RADIUS 

DO 574 I • 1# I A Vc 

CALL MOV A BS ( 930# IY ) 

CALL AOUTS T ( 4# II(I) ) 

CALL AOUTS T ( 5# TA BL E2 ( I ) ) 

IY - IY - 15 

574 CONTINUE 

CALL MOVABS ( 86C# 65C ) 

CALL AOUTSK 5# TABLE2(i) ) 

CALL MOVABS ( 6 GO# 40 ) 

CALL AOUTSK 5# TABLE2(IAVE) ) 

CALL C HRS 1 2 4 2 ) 

CALL MOVABS! 840# 360 ) 

CALL AOUTSK 1# *R • ) 

CALL MOVABS! 38C# 260 ) 

CALL AOUTSK 5# 'ANGLE' ) 

CALL MOVABS! iOO# 5 ) 

IF ( IOP .NE. 11 ) GO TO 580 
CALL aOUTS T ( 1# • I • ) 

CALL CHRS I Z ! 4 ) 

CALL AOUTS T ( 3# 'MAX' ) 

CALL C HRS I Z ! 2 ) 

CALL AOUTS T ( 15# •# R# A DISPLAY' ) 

GO TO 590 

C DISPLAY LABEL 

59C IF ! IOP .NE. 12 ) GO TO 581 

CALL AOUTSK 27# 'I# R# A DISPLAY AT VELOCITY 


5# ANGAVU) ) 

5# ANGAV(KNTll) 

5# ANGAV(KNTl) ) 
5# ANGAVU) ) 


CELL ' ) 
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CALL AOUTSH 3, LCELL > 

GO TO 590 

531 IF ( IOP .NE. 21 > GO TO 582 

CALL AOUTS T ( 32/ 'V, R, A DISPLAY (VELOCITY AT I* » 

CALL CHRSIZ(4> 

CALL AOUTS T ( 3/ 'MAX* ) 

CALL CHRS l Z ( 2 ) 

CALL ACUTS T ( 1/ ')• ) 

GO TO 590 

582 IF ( IOP .NE. 22 ) GO TO 583 

CALL AOUTS T ( i, 'V* ) 

CALL C HRS I Z ( 4 ) 

CALL A OUTS T ( 3/ *M AX' ) 

CALL CHRS I Z ( 2 ) 

CALL AOUTS T ( 15/ ', R, A DISPLAY* I 
CALL CHRS I Z ( 4 ) 

CALL AOUTST< 26/ * AMPLITUDE THRESHOLD AT' ) 

CALL CHRS I Z ( 3 ) 

CALL AOUTSH 4, • I » ' ) 

WRITE ( 2/ 904 ) IVAMP 
REWIND 2 

READ ( 2/ 93i ) IVAMPL 

REWIND 2 

CALL A OUT S T ( 4/ IVAMPL ) 

GO TO 59G 

583 IF ( IOP .NE. 31 ) GO TO 590 

CALL AOUTS T ( 5/ *V < I* ) 

CALL CHRS I Z 1 4 ) 

CALL AOUTS T ( 3/ 'MAX* ) 

CALL CHRS I Z ( 2 ) 

CALL A OUTS T ( 15/ ', R/ A DISPLAY' ) 

590 CONTINUE 

C THREE-DIMENSIONAL DISPLAY ACTIVATED 

CALL THREED ( Rz THETA/ PHI/ ZMIN/ Z MAX/ NPX/ NPY/ NSC ALE/ 

1 NH IDE / NBORD/ NGRID/ Z ) 

CALL A I NST I 3/ ICRIT ) 

CALL ERASE 

CALL MOVA BS ( 100/ 280 I 
595 CALL AOUTST ( 13/ 'ENTER OPTION* ) 

I F( IDAY.GE.238 ) GO TO 592 
C ENTER OPTION FOR UNTRANSLATED DATA 

CALL A I N$ T ( 2/ IOPTN ) 

CALL CORE ( IOPTN/ 2 ) 

READ ( 99 , 602 ) IOPTN 
I OP* I OPTN 

C IF OPTION IS 55 CALL "RAY" 

IF ( IOP .EQ. 55 ) GO TO 71C 
GO TO 596 

C ENTER OPTION FOR TRANSLATED OATA 

592 CALL A INST ( 3/ I OPTN) 

CALL CORE! IOPTN, 3) 

R EAD( 99, 903 ) IOPTN 
IOP-IOPTN 

C IF OPTION IS 555 OR 55 CALL "RAY" 
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IF ( IOPTN .EQ. 555 ) GO TO 710 
5*6 IF( IOPTN.GT. IOC) 10P-I0PTN/10 
IF ( IOP .EO. 55 ) GO TO 710 
I F( IOP.NE.22 )GC TO 597 
CALI MOVABSI IOC, 260 I 

CALL AOUTST ( 43# 'ENTER AMP FOR MAX Vc L OC IT Y . .. ( 3 OIGITS)' ) 
C ENTER AMPLITUDE THRESHOLD 

CALL AI>N$T{ 3, I V AMP ) 

CALL CORE ( I VAMP, 3 > 

READ ( 99, 903 ) IVAMP 
C READ RAW DATA 

READ(3,9b4, END-587 > I FRM, XY( 4 ) , I AMP 
XW< 1>*1 

C ASSIGN MAXIMUM VELOCITY 

00 585 1-2,100 
IF( IAMP(I).GT. IVAMP) XW<i>«I 
585 CONTINUE 

DO 587 K-2, NN 

RE AD 13, 964/ END -5 87) IFRM,XY(4),IAMP 
XW< I )*1 

DO 586 J-2, 100 

IF( lAMPIJI.GT. IVAMP) XW<I)«J 

536 CONTINUE 

537 CONTINUE 
REWIND 3 

597 CONTINUE 

C IF OPTION IS 99 OR 999 RETURN TO CALLING PROGRAM 

IF ( IOPTN .EQ. 99 ) RETURN 
I F ( I0PTN.EQ.99t) RETURN 
IF< IOPTN. EQ. 999) RETURN 
CALL ERASE 
GO TO 548 
710 CONTINUE 

CALL RAY( NFRM1, ANG 1 , NFRM2, ANG2, NS CAN, ISCAN1, IOAY, IRUN, 

1 NT APE , ITAPES, IC ELL , I SC AN, ITMIN, ITMAX ) 

GO TO 1 


901 

FORMAT 

( 

11 

) 


902 

fdrmat 

( 

12 

) 


903 

FORMAT 

( 

13 

) 


90 4 

FORMAT 

( 

14 

) 


90 5 

format 

( 

15 

) 


906 

format 

( 

16 

) 


90 8 

FORMAT 

( 

18 

) 


921 

FORMAT 

( 

F 5. 

2 

) 

922 

format 

< 

F 4 . 

2 

) 

931 

FORMAT 

( 

A4 

) 


932 

format 

I 

2A4 

) 


952 

FORMAT 

( 

F 6 . 

2 

) 

95 i 

FORMAT 

1 

20 A 4 

) 

999 

RETURN 






END 
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****** ************* ********* ************ ** ************** ** ****** ***** 

* 

SUBROUTINE RAY * 

* 

THIS SUBROUTINE DISPLAYS ONE "RAY** OF DATA IN THREE ♦ 

DIMENSIONAL FORM * 

* 

********** *********************************************************** 


SUBROUTINE RAY ( NFRMi* ANG1* NFRM2* ANG2* NSC A N* ISCAN1* IDAY, 
1 iRUN* NTAPE* ITAPES* ICELL* ISCAN* ITM IN* ITMAX ) 

DIMENSION I AM P ( 1 04 ) * XY14I* INI2) 

DIMENSION ZtiCC*iOCI* NGRID14J 
DIMENSION TITLE120) 

DIMENSION NFRM 1 1 60 ) * ANGildd* NFRM2I6GI* ANG2I6C')* ITAPESI2J 

DQUBLEPREC IS ION LFRM1, AANG1* LFRM2* AANG2* LTAPEi * LTAPE2 

DOUBLEPRECISION ANGAV I 50 ) * TABLE2150I 

DOUBLEPRECISION ZAVE2* ZHIGH2 

COMMON / HARP / ICAL*NSCOPE 

CALL CHRSIZI 3 ) 

NTAPE3 • 0 
ISAME * C 
IRET « 0 
NPX ■ ICO 

C SEE WHICH TAPE NEEDS TO BE READ 

IF ( ISCAN .GE. ISCANI ) NTAPE*3 
NSCOPE * 3 

C FIND SCAN L IMITS 

WRITE ( 2* 904 ) ISCAN 
REWIND 2 

READ ( 2* 931 ) LSCAN 
REWIND Z 
CALL ERASE 

IREAD«i 

CALL CHRSIZI3) 

C DISPLAY TO USER WHICH SCAN HE IS IN 

CALL MQVA85 ( 100* 750 ) 

CALL AOUTSTI 6 * * SCAN #• ) 

CALL AOUTSTI 4* LSCAN > 

CALL MOVABSI IOC* 720 ) 

CALL AQUTSTI17 * ’ENTER ANGLE ( R EAL ) ' ) 

C ENTER ANGLE SIZE 

CALL AINSTI 4 * A NG ) 

CALL CORE! ANG* 4 ) 

READ ( 99, 921 ) ANG 

921 FORMAT I F4.2 ) 

CALL MOVABSI 100* 690 ) 

CALL AOUTSTI 18* ’ENTER SPRE ADI RE AL I ' I 
C ENTER SPREAD 

CALL AINSTI 4* SPR ) 

CALL CORE I SPR* 4 ) 

READ I 99, 922 ) SPR 

922 FORMAT I F4.3 ) 

C IF ANGLE IS ZERO USE ANGLE RANGE FROM CALLING PROGRAM 



0& PA GE IS 

eoor qcaijjy 


c 


202 

c 

911 
210 

912 

215 


C 


22 0 
C 

299 


400 


40 4 


402 


IF 

C 

ANG .EO. 0.0 

) 

I ANG « ITMIN 

IF 

( 

ANG .EQ. 0.0 

) 

JANG * ITMAX 

IF 

( 

ANG .EQ. C.O 

) 

GO TO 202 


ASSIGN ANGLE LIMITS 
I ANG * ANG - SPR 
JANG « ANG ♦ SPR 
SANG * 90.0 
CALL ERASE 

CONVERT NUMERICAL DATA TO CHARACTER FORM 
READ { i, 911 ) IDUH/ JOUM/ KFRM, MFRM 

FORMAT ( 12, II# 16/ 16 > 

LFRM ■ KFRM 
KFRM - LFRM 

READ ( 1/ 912 ) I SC N/ IDRC/ TAVE/ LFRM 

FORMAT ( 12/ 11/ F 6. 2/ 16 ) 

J SC N • ISCN 

IF ( I SC N .Nfc. ISCAN > GO TO 210 
IAVE - TAVE 

IF ( IAVE .Gt. I ANG .AND. IAVE .LE. JANG ) GO TO 220 
SANG » TANG 
KFRM • LFRM 

READ ( 1/ 912 > ISCNz IDRC/ TAVE/ LFRM 

IS ANGLE GIVEN 8Y USE IN SCAN RANGE? 

IF ( JSCN .£0. I SC N ) GO TO 215 

CALL AOUTS T ( 2 7 / 'ANGLE NOT IN SCANNING RANGE' ) 
CALL AINST( 1/ DUMMY ) 

REWIND 1 
RETURN 
CONTINUE 
IRET « 0 

READ DATA 


IF 

( 

ntape 

.EQ. 3 

) 

CALL 

GET2 ( I FRM/ 

XY, 

ITM/ 

I AMP/ 

IRET 

) 

IF 

( 

NTAPE 

• LE. 2 

) 

CALL 

GET ( IFRM/ 

XY/ 

I TMz 

I AMP/ 

IRET 

I 

IF 

( 

I FRM . 

LT. KFRM 

) GO 

TO 299 






IF 

( 

I FRM . 

EQ. KFRM 

) GO 

TO 400 






IRET 

* 1 










IF 

( 

NTAPE 

-LE . 2 

) 

CALL 

GET (IFRM/ 

XY/ 

ITM/ 

I AMP/ 

IRET 

) 

IF 

I 

NTAPE 

.EQ. 3 

) 

CALL 

GET2I IFRM/ 

XY, 

ITM/ 

IAMP/ 

IRET 

) 


IRET - 0 
GO TQ 299 

IF l IDRC .EQ. 0 I GL TO 412 
IF ( IDRC .EQ. 1 ) GL TO *02 

CALL AOUTSTI 28/'ERROR IN DIRECTION INDICATOR' ) 

CALL AlNSTI 1/ DUMMY ) 

CALL ERASE 
GO TO 999 

READ ( 1/ 912/ E NO = 999 ) ISCN/ IDRC/ TAVE/ LFRM 
IF ( ISCN .NE. ISCAN ) GO TO 999 
IAVE » TAVE 

IF ( IAVE .GE. I ANG .AND. IAVE .LE. JANG ) GO TO 402 

GO TO 999 

NPY » LFRM - KFRM 

K » NPY 

IMAXZ - 0 
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C 


C 


406 

407 


C 

C 


C 


c 

c 


KQUNT * K3UNT ♦ 1 
READ DATA 
00 407 I » 1, NPY 

IF ( NTAPE .LE. 2 ) CALL GET (IFRM# XY# UK, IaMP# IRET > 

IF t NT APE .EQ. 3 ) CALL GET2IIFRM# XY, ITM# IAMP# IRET ) 

I F t t.EQ.i) A»XYI 3> 

AA»XY( 3 ) 

FIND MAXIMUM AMPLITUDE 
DO 406 J * 1# KPX 

IF ( IMAXZ .LT. UMPIJ) ) IMAXZ • IAMPUJ 
HK,Ji » IAMPIJ+2) 

K • K - 1 
CONTINUE 
CALL ERASE 

CHANGE ANGLE TO TRUE ANGLE 
AANG1 - ANG1IISCAN) * 2.0 
AANG2 * ANG2IISCAN) * 2.C 

CONVERT NUMERICAL OATA TO CHARACTER FORM 
WRITE ( 2* 962 ) NFRM1 ( I SCAN )# AANG1# NFRM2I I SCAN) # 

1 AANG2 # I DA Y# I RUN# ITAPES# ISCAN, ICELL 

REWIND 2 

REAO ( 2# 963 ) LFRMi# A ANG1 * LFRM2# AANG2# LOA Y» L RUN# LTAPE 1 # 

1 LTAPE2# LSCAN# LCELL 

REWIND 2 

LABEL DISPLAY 
CALL C HRS 12(2) 

CALL MOVABSI 100* 750 ) 

CALL AOUTSTI 32# ' NASA PROJECT DUST DEVIL DAY • I 

CALL AOUTSTI 3# LOAY ) 

CALL AOUTST ( 10# ' RUN « • ) 

CALL AOUTSTI 2# L RUN > 

IF INTAPE. EQ. 3) CALL AOUTST 1 16# * TAPE 2 OF 2-*) 

IF INTAPE. EQ. 2) CALL AOUTST 116#' TAPE 1 OF 2-') 

IF I NTAPE .EQ. i ) CALL AOUTSTI 13# • TAPE # 1-' ) 

IF INTAPE . E Q. 3 ) CALL AOUTST (5#LTAPE2) 

IF INTAPE. LE. 2) CALL aOUTST I5#LTAP6l» 

CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 
CALL 

CALL 


CHRSIZI3) 

MOVABS I 80# 725 I 
AOUTSTI 7# 'SCAN * • ) 

AOUTSTI 2# LSCAN ) 

AOUTSTI 13# * FRAME # * ) 

AOUTSTI 6# LFRM1 ) 

AOUTSTI 9# • ANGLE • ) 

AOUTSTI 6# AANG1 ) 

AOUTSTI 18# ' TO FRAME « ' ) 

AOUTSTI 6# LFRM2 ) 

AOUTSTI 9# • ANGLE ' I 

AOUTSTI 6# AANG2 > 

AOUTSTI 25# • CRITICAL AMP AT CELL • > 

AOUTST I 3# LCELL ) 

MOVABSI 70# 50C > 

MOVABSI 10G# 3C0 ) USE THIS TO SEE GRAPH FROM 

DIFFERENT ANGLE 

AOUTSTI 8# * RADIUS' ) 
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C 


C 


C 


C 

C 


C 


C 


CONVERT NUMERICAL DATA TO CHARACTER FORM 
WRITE ( 2# 959 ) A/ AA 
REWIND 2 

READ C 2# 939 ) A, AA 

REWIND 2 
CALL CHRS 12(41 
CALL MOVABS ( 0# 340) 

CALL MOVABS ( 38G# 140 ) 

CALL A0UTST!4#AI 
CALL MOVABS (360> 64C ) 

CALL M0VABS<0# 450 ) 

CALL AOUTST ( 4# AA ) 

CALL CHRS 12(3) 

LABEL FOR GRAPH 
CALL MOVABS ( 30C.» 140 ) 

CALL MOVABS ( 75C# 200) 

CALL AOUTST! 22# * VE LOC I TY . . AMP< 3 - 1021* ) 
CALL MOVABS! 700, 700 ) 

CALL AOUTST ( 7, ’FRAME A ' ) 

WRITE! 1, 908 ) KFRM 
REWIND 2 

READ ( 1, 932 ) IN 
REWIND 2 

LABEL FOR GRAHP 
CALL AOUTST! 8, IN > 

CALL AOUTST! 4, ' TO’ ) 

WRITE ! 2# 908 ) LFRM 
REWIND 2 

READ ( 1, 932 ) IN 
REWIND 2 

CALL AOUTST ( 8, IN ) 

CALL MOVABS! 700# 670 ) 

CALL AOUTST ( 7# ’ANGLE »’ ) 

WRITE ! 1, 95 1 ) TAVE 
REWIND 2 

READ ( 2# 932 ) IN 

REWIND 2 

CALL AOUTST ( 8# IN ) 

CALL MOVABS! 7C0# 640 I 

CALL AOUTST! 19# ’MAXIMUM AMP VALUE » • ) 
WRITE 1 2# 904 ) IMAXZ 

REWIND 2 

READ I 1, 931 ) I MAXZZ 

REWIND 2 

CALL AOUTST! 4# IMAXZZ ) 

SET UP THE OPTION FOR 3-0 DISPLAY 
R • 500.0 
ICAL • 0 
THETA ■ 30.0 
THETA --30.0 
PHI • 30.0 
ZMIN • 0.0 
ZMA X » 20.0 
NPX ■ 100 


USE THIS FOR DIFFERENT ANGLE 
USE THIS FOR DIFFERENT ANGLE 

USE THIS FOR DIFFERENT ANGLE 


C 


USE THIS FOR DIFFERENT ANGLE 



NHIOE « 1 
NSCALE « 2 
NBORD * 0 
NGRID(i) » C 
NGR ID ( 2 ) » 0 
NGRIDC3) • 0 
NGR ID I 4 ) • 0 
ZMN » Z<1»1) 

ZMX - Z(l*i> 

C FIND MAXIMUM ANO MINIMUM DATA 

C FOR SCALING FACTORS 

DO 451 I * 1* NPY 
DO 451 J * 1* NPX 
IF ( ZMN .GT. Z<I*J> ) ZMN « Z(l*J> 

IF ( ZMX .LT. Z(I*J» ) ZMX » ZU«J) 

451 CONTINUE 

C SCALE DATA 

DO 452 I • 1 * SPY 
DO 452 J * 1 * NPX 

Z(I*J) « ( Z(I*J> - ZMN » * ( 10Q.0 / ( ZMX - ZMN I ) 

452 CONTINUE 

C PLOT 3-D GRAHP 

CALL THREEDI R* THETA* PHI * ZMIN* ZMAX* NPX* NPY* NSCALE* 

1 NHIDE* NBORD* NGRIO* Z > 

KFRM » LFRM 


C CONVERT NUMERICAL DATA TO CHARACTER FORM 


414 READ 

( 1* 912* 

END-999 

) 

ISCN* IDRC* 

TAVE* 

LFRM 

IF < 

ISCN .NE. 

ISCAN ) 

GO 

TO 999 



I AVE 

■ TAVE 






IF ( 

IAVE .GE. 

I AN G .AND. 

IAVE .LE. 

JANG ) 

GO TO 412 


GO TO 999 

412 NPY » LFRM - KFRM 
IMAXZ * 0 

KOUNT * KOUNT ♦ 1 
C READ DATA 

DO 415 I » 1* NPY 

IF ( NTAPE .EQ. 3 > CALL GET2UFRM* XY* ITM* IAMP* IRET ) 

IF I NTAPE .LE. 2 ) CALL GET I I FRM* XY* ITM* IAMP* IRET ) 
IFft.EQ.l) A«XY(3> 

AA«XY<3> 

DO 413 J « 1* NPX 

IF ( IMAXZ .LT. IAMP(J) ) IMAXZ « IAMPIJJ 

413 Z(I*J) • IAMPIJ+2) 

415 CONTINUE 

CALL ERASE 

AANG1 ■ ANGKISCAN) * 2.C 
AANG2 « ANG2 ( I SC AN ) * 2.0 

WRITE < 2* 962 » NFRM 1 ( ISCAN)* A A N G 1 * NFRM2I I SCAN) * 

1 AANG2* IDAY* IRUN* ITAPES* ISCAN* ICELL 

962 FORMAT ( 16* F6.3* 16* F 6.3* 13* 12* 215* 12* 13 ) 

REWIND £ 

READ ( 2* 963 ) LFRM1* AANG1* LFRM2 * A4NG2 * LOAY, LRUN* LTAPE1* 
I LTAPE2* LSCAN* LCELL 

963 FORMAT ( 4A6* A3* A2* 2A5* A 2* A3 ) 
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C 


C 


C 


REWIND 2 

CALL CHRSIZI2) 

CALL MOVABS! 10J* 750 ) 

CALL AOUTSTI 3 2* 'NASA PROJECT OUST DEVIL DAY * ) 

CALL AOUTSTI 3* LDAY > 

CALL AOUTST ( 1C* • RUN # • ) 

CALL AOUTSTI 2* LRUN ) 

IF INTAPE. EQ. 3) CALL A0UTST116*’ TAPE 2 OF 2-*) 

IF INTAPt.EQ.2) CALL AOUTST 116*' TAPE i OF 2-') 

IF C NT APE .£0. 1 ) CALL AOUTSTI 13* • TAPE # l-’ ) 

IF INTAPE . EQ. 3 ) CALL AOUTST 15*LTAPE2I 
IF I NT APE .L E . 2 ) CALL AOUTST 15*LTAPE1) 

CALL CHRSIZ<3> 

CALL MOVABS ( 80* 725 ) 

CALL AOUTSTI 7* ’SCAN « • ) 

CALL AOUTSTI 2* LSCAN ) 

CALL AOUTSTI 13* * FRAME 8 • ) 

CALL AOUTSTI 6* LFRM1 ) 

CALL AOUTSTI 9* • ANGLE • ) 

CALL AOUTSTI 6* AANG1 ) 

CALL AOUTSTI 18* • TO FRAME # • ) 

CALL AOUTSTI 6* LFRM2 ) 

CALL AOUTSTI 9* • ANGLE * ) 

CALL AOUTSTI 6* AANG2 ) 

CALL AOUTSTI 25* * CRITICAL AMP AT CELL * > 

CALL AOUTST I 3* LCELL ) 

CALL MOVABS I 70* 500 ) 

CALL MOVABS! ICO* 300 ) USE THIS FOR DIFFERENT ANGLE 

CALL AOUTST! 8* * RADIUS’ ) 

WRITE I 2* 959 ) A* AA 

REWIND 2 

READ I 2* 939 ) A. AA 

REWIND 2 

CALL CHRSIZI4) 

CALL MOVABS I 36 C* 640 ) 

CALL MOVABS! 0* 450 ) USE THIS FOR DIFFERENT ANGLE 

CALL AOUTST I 4* A ) 

CALL MOVABS I 0* 340 ) 

CALL MOVABS! 38C* 140 ) USE THIS FOR DIFFERENT ANGLE 

CALL A0UTST(4*AA> 

CALL CHRSIZI3) 

CALL MOVABS! 30C* i40 > 

CALL MOVABS! 750* 200) USE THIS FOR DIFFERENT ANGLE 

CALL AOUTSTI 22* ' VE LOC IT Y . . AMP I 3 - 1C2)» ) 

CALL MOVABS! 70C* 700 ) 

CALL AOUTST I 7* ’FRAME #* ) 

WRITE! 2* 908 ) KFRM 
REWIND 2 

READ I 2* 932 > IN 

REWIND 2 

CALL AOUTSTI 8* IN ) 

CALL AOUTSTI 4* • TO* ) 

WRITE I 2* 908 ) LFRM 
REWIND 2 



61 


C 


C 


461 

C 


462 

C 


99 9 
90 1 


READ I 2# 932 I IN 
REWIND 2 

CALL AOUTST ( 8 » IN ) 

CALL MOV A 8S ( 7QC# 670 ) 

CALL AOUTST ( 7, ‘ANGLE «• ) 

WRITE ( 2# 952 I TAVE 
REWIND 2 

READ ( 2, 932 ) IN 

REWIND 2 

CALL AOUTST ( 6 # IN ) 

CALL M0VA8SI 70C# 640 ) 

CALL AOUTSTl 19, 'MAXIMUM AMP VALUE 
WRITE < 2# 904 > IMAXZ 
REWIND 2 

READ ( 2# 93 i ) IMAXZZ 
REWIND 2 

CALL AOUTST ( 4# IMAXZZ ) 

OPTIONS FOR 3-D DISPLAY 
NBORD » 0 
NGR ID 11 > « 0 
NGR I D ( 2 ) ■ C 
NGR ID ( 3 ) • 0 
NGR ID ( 4 I » C 
R * 500,0 
ICAL • 0 
THETA » 30.0 
THETA *-30. 0 
PHI • 30.0 
ZMIN « 0.0 
ZMAX « 20.0 


°*^0OR of AOE K 
°° R QUALITY 


USE THIS FOR DIFFERNT ANGLE 


NPX » ICO 
NHIDE • 1 
NSCALE » 2 
ZMN » Z<1#1> 

ZMX • Z(i#l> 

FIND SCALING FACTOR 
DO 461 I » 1# NPY 
DO 461 J ■ l, NPX 

IF I ZMN .GT. Z(I#J> ) ZMN » ZU#J> 

IF ( ZMX .LT. Z(I#J> ) ZMX « Z(I#J> 

continue 

SCALE DATA 
DO 462 1*1# NPY 

DO 462 J * 1 # NPX 

* < Z(I#J) - ZMN J * ( 100.0 / ( ZMX - ZMN I > 

CONTINUE 

PLOT 3-D GRAHP 

CALL THREE D ( k # THETA# PHI# ZMIN# ZMAX# NPX# NPY# NSCALE# 

1 NHIDE# NBORD# NGR ID# Z > 

KFRM » LFRM 


GO TO 404 


REWIND 1 
RETURN 

FORMAT < Ii ) 
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FORMAT 

( 

13 ) 

904 

FORMAT 

< 

14 > 

905 

FORMAT 

< 

15 ) 

906 

FORMAT 

< 

16 1 

906 

FORMAT 

( 

18 ) 

931 

FORMAT 

( 

A 4 ) 

932 

FORMAT 

< 

2A4 ) 

939 

FORMAT 

( 

2A4 ) 

952 

FORMAT 

( 

F 6 . 2 

951 

FORMAT 

( 

2CA4 

959 

FORMAT 

END 

1 

2F4.0 
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SUBROUTINE T HR EED 

PLOTS THREE-D IMENS IONAL FIGURE USING LIBRARY PROGRAMS 


• * 

* 

* 

* 


SUBROUTINE THREtOI R » THETA, PHI, ZMIN, l MAX# NPX# NP Y# 

I NSCALE# NH IDE# NBORD# NGRID# Zi ) 

R-- IS THE DISTANCE FROM THE CENTER OF THE PLOT WHICH YOU WISH 
TO PUT THE CAMERA 

THETA-- IS THE ROTATION IN DEGREES IN THE XY PLANE WHICH YOU 
WISH TO PUT THE CAMERA 

PHI-- IS THE ANGLE OF ELEVATION OF THE CAMERA FROM THE XY PLANE 
ZMIN AND ZM AX ARE THE LOWER AND UPPER LIMITS OF THE Z VALUES 
NPX IS NUMBER OF ARRAY POINTS IN X DIRECTION 
NPY IS NUMBER OF ARRAY POINTS IN Y DIRECTION 
NSCALE DEFINES THE SCALING METHOD FOR SCALING Z ARRAY TO 
BOX HEIGHT. NSCALE MAY TAKE ON THE FOLLOWING VALUES: 

1 - THE RANGE OF Z FROM LOWER TO UPPER CUTOFF GIVEN IN 

ZMIN AND ZMAX EXACTLY FITS BETWEEN TOP AND BOTTOM HALF 
OF BOX. 

2 - SCALE TO FIT ENTIRE RANGE OF Z ARRAY INTO BOX. 

3 - SCALE SO THAT Z*0# KEEPS SURFACE WITHIN BOX. 

4 - SAME AS l BUT PUT CUTOFF VALUES ON BOUNOARY 

5 - SAME AS 3 BUT PUT CUTOFF VALUES ON BOUNDARY. 

6 - USE CUTOFF VALUES FOR SCALING 
NH IDE - IF ZERO DRAW ALL POINTS# IF NONZERO DO NOT DRAW 

HIDDEN POINTS. 

•NGRID - IF NGPIO(i) • 0# NO COORDINATE GRIDS ARE DRAWN. 

OTHERWISE# NCR I D ( 2 ) * NUMBER OF LINES FOR X AXIS 

NGR ID ( 3 ) * NUMBER OF LINES FOR Y AXIS 

NGRIGI4) * NUMBER OF LINES FOR Z AXIZ 

COMMON / HARP / ICAL#N$COPE 

DIMENSION I BUF ( 7000 )# Z 1( IOC# 100 ) # I FORM 1 (10) 

DIMENSION m0C#i00l#XYLIM(2#6l#CAMLQC( 32 ) # NGR ID ( 4 ) # 

1 LAB(2C># IFORM(IG) 

DATA I FORM / ' ( 20 • # • F 5 .0 • # • ) •#• •#« •#• •#• •# 

1 * •# • '# • • / 

NSC OP E * 0 
ICOUNT«C 
N S K I P * C 
K T R A C K * C 

I IFIKTRACK.NE.O ) CALL TR I P I N ( NSC ALE # Z MI N # Z M A X# NH 10E # NBORD# NGR I D# K TR 
#ACK#KLOG#XCAM# YCAM# ZCAM# NPLOTS# Z# Zl#NPXi#NPYi»ICAL) 

I F ( ICAL.EO.i )IC0UNT*ICGUNT+1 

IFIKTRACK.EQ.999) GO TO 41 

IF ( KTRACK .EQ. 0 > GO TO 90 

XMIN-XMINi 

XMAX-XMAXi 

YM IN* YMI N1 



<~> r> o o r> o 


YMAX- YMAX1 
NPX-NPX1 
NPY-NPYi 
90 CONTINUE 

IFIKTRACK.NE.O ) GO TO 33 
XCAM « R. 

YC AM * THETA 
ZCAM « PHI 
00 92 K=1,NPY 
DO 92 J*1,NPX 
92 Z C K> JI-Z1I K, J) 

33 IFIKLOG.EO.i > GU TO 96 
XMIN1-XMIN 
XMAX1-XMAX 
YMIN1 -YMIN 
YMA XI • YMAX 
NPX1-NPX 
NPYl-NPY 
96 CONTINUE 

HEIGHT - 7.50 
WIDTH - 10.50 




MINIMA AND MAXIMA NOT PROVIDED 

BOT-ZI i/i > 

TOP-BOT 

DO 30 I * i > NP Y 
DO 30 J « l > NP X 

I F I SOT .GT. Z(I*J>> BOT » Z(I»J) 
I F I TOP .IT. ZII,J)) TOP - ZII/J) 
30 CONTINUE 
XMIN * BOT 
XMAX » TOP 
YMIN « BOT 
YMAX * TOP 

32 IFIZMIN .NE. ZMAX! GC TO 29 
ZMIN « BOT 
ZMAX • TOP 
29 CONTINUE 


NOW STORE THESE VALUES AS NEEDED FOR FUTURE USE 


XYl IM( 1# 1 1 
X YL I M I 2> 1 ) 
XYLIMI 1*2) 
XYLIMC2>2) 
XYLIM(i>3) 
XYLIMI 2,3> 
X YL I M I i » 4 ) 
XYLIMI 2/ 4) 
XYLIMU/ 5) 
XYLIM(2»5> 
XYLlM(l*6) 
XYLIMI2/6) 


XMIN 

XMAX 

YMIN 

YMAX 

ZMIN 

ZMAX 

ZMIN 

ZMAX 

W IDTH 

HEIGHT 

0.5 

0.5 
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X VALUES 
- % 14 I 


*, 14 , /, 


INCLUDING HID06N POINTS - DR4WN' ) 


NOT DRAWN' ) 


WR I Tt ( 6, 7 ) NPX,NPY 
7 FORMAT!/,' NUMBER OF 
1' NUMBER OF Y VALUES - 
IF1NHIDE ) 6/8, 10 
WRITE 16,9) 

FORMAT!' ALL PCINTS - 
GO TO 12 
WRITE 16, *1 ) 

FORMAT!' HIDDEN POINTS 
IF1NB0RD) 13,15,13 
WRITE lb, 14) 

FORMAT!' BORDER WILL BE 
GO TO 17 
WRITE! 6, lb) 

FORMAT!' NO BORDER AROUND 
I F ! NGR ID ( l ) ) i8,ld,2C 
WRITE !b, 19) 

FORMAT!' GRID LINES WILL NOT 
GO TO 22 

WRITE ! 6, 21) NGR I 012 ) , NGR I D I 3 ) , NGR I D 1 4 I 
FORMAT!' GRID LINES SELECTED 
15X, 14,' GRID LINES WILL BE DRAWN ALONG 

25X, 14,' GRID LINES WILL BE DRAWN ALONG 

35 X, Ia, ' GRID LINES WILL BE DRAWN ALONG 

22 WRITE!b,24) 

DO 25 I = 1,NPY 

24 FORMAT!///, • ///// INPUT MATRIX FOLLOWS 

25 WRITE 16, IFORM) ( Z( I , J ), J-l, NPX) 
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DEFINE POSITIONS 


DEFAULT OPTIONS REMOVED 

CHANGING CAMERA VALUES TO BE ENTERED IN TERMS OF R, THETA, PHI 
PI * 3.141 5926 
YCAM* PI*YC AM/180. 

ZCAM* P I *ZC AM / 1 8G. 

THETA* YCAM 

YCAM* XC AM ♦ SIN (YCAM) 

ZCAM* XC AM ♦ SIN(ZCAK) 

XCAM * XC AM * CGS(THcTA) 

C AMLOC! 1 ) « XCAM 
C AMLOC ( 2 ) « YCAM 
C AMLOC ( 3 ) « Z CAM 
DO 27 I * 4,32 
27 C AMLOC! I ) ■ C.O 
CAMERA AIM 

C AML OC ( 9 ) * ( XMAX-XMINW2. 

CAMLOC(iO)* (YMAX- YMINI/2. 

C AMLOC ( li ) * 0.0 

READY TO PLOT 


CALL CHRSIZI4I 

CALL VWIND0(G. ,11.,C.,8.5) 
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CALL PURJOYI Z, NPX,NPY,iOC*CAMLOC>XYL W>Oj NS CAL E> NHI DE# N80RD# NGR I D) 

LOOP FOR NEXT DATA SET 

IF( NSCOPE.EQ.O ) KTRACK* KTRACK +1 
I F ( NSCOPE.NE .0 ) KTRAC K*0 
GO TO i 

41 IF(NPL0TS.EQ.999) GO TO 28 
28 CONTINUE 
RETURN 
END 
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